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Summary  of  Achievements 

In  the  three  years  of  this  grant,  we  have  made  significant  progress  towards  the  goal  realizing  scalable  quantum  bits  for 
information  processing,  using  charge  states  in  superconducting  electronics.  Specifically,  we  have  made  a  single  qubit,  based  on  the 
single  Cooper-pair  box  (CPB),  which  has  a  phase  coherence  quality  factor  greater  than  100,  and  a  very  long  inelastic  lifetime 
corresponding  to  an  intrinsic  quality  factor  greater  than  100,000.  These  results  were  obtained  by  performing  the  first  ever  charge 
measurements  using  a  radio-frequency  single-electron  transistor  (RF-SET),  integrating  it  with  the  charge  qubit,  and  performing  both 
precision  microwave  spectroscopy  and  sub-microsecond  resolved  readouts  of  the  qubit.  Our  results  indicate  that  single-shot 
measurements  of  individual  qubits  with  this  technique  will  be  possible.  These  results  have  been  accepted  for  publication  in  Physical 
Review  Letters.  Moreover,  we  have  studied  the  backaction  of  this  type  of  charge  measurement  with  an  SET,  both  experimentally  and 
theoretically,  which  is  relevant  to  several  other  schemes  for  quantum  computation  that  invoke  similar  measurement  strategies. 
Significant  effort  was  made  in  building  fabrication  expertise,  experimental  techniques,  and  apparatus.  This  project  is  continuing  under  a 
new  three-year  grant,  during  which  we  use  these  capabilities  to  further  investigate  the  mechanisms  and  limits  of  decoherence  in  the 
CPB  qubits,  test  the  models  of  quantum  measurement  by  the  RF-SET,  and  investigate  the  entanglement  of  two  CPB  qubits  using 
techniques  analogous  to  cavity-QED  in  atomic  physics. 

•  Developed  fabrication  process  at  Yale  for  Al-based  SETs  using  electron-beam  lithography. 

•  Began  finite-element  full-wave  and  electrostatic  simulations  of  actual  chip  geometries. 

•  Revised  chip  designs  for  better  control  of  electromagnetic  coupling  to  gates. 

•  Measured  equilibrium  charge  state  of  Cooper-pair  box  and  observed  2e-periodic  Coulomb  staircase. 

•  Performed  CW  spectroscopy  of  charge  superpositions  in  box  with  microwave  signal  on  gate. 

•  Tuned  parameters  of  the  CPB  qubit  using  applied  gate  voltages  and  magnetic  fields. 

•  Accurately  determined  all  parameters  of  the  qubit  Hamiltonian  by  spectroscopy. 

•  Determined  ensemble-averaged  decoherence  time  (T2*  ~  1  ns)  from  width  of  photon  resonances. 

•  Observed  dependence  of  CPB  ground-state  polarization  on  bias  conditions  of  readout  SET. 

•  Using  time-gated  RF-SET  measurement,  observed  mixing  of  states  by  measurement  backaction. 

•  Measured  long  inelastic  lifetime  of  CPB  excited  state  (Ti=1.3  gs),  under  conditions  of  continuous  measurement. 

•  Separated  effects  of  electric  and  magnetic  noise  on  relaxation. 

•  Estimated  that  inelastic  lifetime  is  sufficiently  long  to  allow  single-shot  quantum  state  readouts  using  the  RF-SET. 

•  A  manuscript  on  the  decoherence  and  lifetime  measurements  have  been  accepted  by  Physical  Review  Letters. 
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Summary  of  Technical  Approach 

The  goal  of  the  project  was  to  investigate  the  limits  on  coherence  and  control  of  quantum  states  in  the  “single  Cooper-pair 
box,”  which  is  a  candidate  system  for  performing  quantum  computation.  These  nanoelectronic  devices  utilize  charging  effects  in  small 
capacitance  Josephson  junctions,  operated  at  millikelvin  temperatures.  Superpositions  of  two  different  charge  states  on  these 
junctions  serve  as  the  qubit,  which  would  form  the  basis  of  a  future  quantum  computer.  Experiments  in  Japan  by  Nakamura  et  al. 
have  shown  that  some  degree  of  coherent  state  control  is  possible,  and  they  achieved  coherent  single-bit  rotations  (Rabi  oscillations) 
between  states.  However,  the  decoherence  processes  and  limits  in  these  devices  are  not  yet  understood  in  detail,  and  these  questions 
were  the  main  focus  of  the  current  project. 

Our  approach  was  to  utilize  the  RF-SET  electrometer,  whose  high  speed  and  good  charge  sensitivity  allow  a  direct 
determination  of  the  island  charge  of  the  Cooper-pair  box.  Since  the  fabrication  used  for  the  RF-SET  is  identical  to  that  of  the  Cooper- 
pair  box,  the  electrometer  can  be  integrated  on-chip  (see  Fig.  1)  to  probe  the  potential  of  the  box  capacitively.  By  fast  control  of  the 
microwave  signal  that  interrogates  the  RF-SET,  the  measurement  can  be  turned  off  during  the  coherent  state  manipulations,  and  then 
the  final  state  can  be  sampled  on  sub-microsecond  timescales.  In  its  quiescent  state,  the  RF-SET  adds  only  a  small  additional 
coupling  to  the  electromagnetic  environment.  With  the  RF-SET  charge  measurements,  we  have  shown  that  it  will  be  possible  to 
approach  the  ideal  limit  of  a  single-shot  measurement  on  the  state  of  a  single  qubit. 


Figure  1 .  A  circuit  diagram  of  a  Cooper-pair  box  qubit,  capacitively  coupled  via  Cm  to  an  RF-SET  electrometer,  showing  the  L- 
C  resonant  tank  circuit. 

These  devices  are  fabricated  using  a  now-standard  fabrication  process,  utilizing  direct-write  electron-beam  lithography  and  a 
multiple  angle  evaporation  of  aluminum  films.  By  oxidizing  between  depositions,  it  is  possible  to  realize  high-quality  AI/AIOx/AI  tunnel 
junctions  with  characteristic  sizes  down  to  about  50  nm.  At  the  millikelvin  temperatures  of  operation,  the  aluminum  is  a 
superconductors,  with  nearly  all  the  conduction  electrons  bound  up  as  Cooper-pairs.  This  allows  one  to  realize  a  solid-state, 
microelectronic  circuit  which  behaves  as  an  effective  spin  V2  particle,  and  thus  a  candidate  qubit.  These  qubits  can  have  the 
parameters  of  their  Hamiltonian  tuned  via  the  application  of  magnetic  and  electric  fields,  and  their  quantum  states  can  be  manipulated 
by  fast  signals  on  a  capacitively-coupled  gate  lead  that  comes  to  room  temperature.  The  samples  for  these  experiments  were 
fabricated  by  our  collaborators  in  Prof.  Per  Delsing’s  group  at  Chalmers  University  in  Sweden,  under  a  subcontract. 

Several  different  types  of  measurements  were  performed  using  the  same  configuration  of  devices.  First,  continuous-wave 
spectroscopy  of  the  box  was  performed,  under  the  conditions  of  a  continuous  measurement  by  the  RF-SET.  From  this  data,  we  could 
extract  the  parameters  of  the  box  Hamiltonian,  verify  the  coherent  two-level  physics  of  our  circuit,  and  obtain  a  worst-case  estimate  of 
the  decoherence  time.  The  intrinsic  high  speed  of  the  RF-SET  then  allowed  a  measurement  of  the  dynamics  on  the  sub-microsecond 
timescale.  By  measuring  the  decay  time  of  an  excited-state  population,  we  could  observe  the  inelastic  relaxation  rate,  which  was  quite 
slow,  and  approached  the  limit  expected  by  spontaneous  emission.  This  was  the  first  measurement  of  this  timescale  in 
superconducting  charge  qubits,  and  very  encouraging  for  future  developments.  Other  dynamical  experiments  were  also  performed, 
with  somewhat  less  success.  The  signal-to-noise  ratio  of  a  single-shot  measurement,  based  on  our  measured  relaxation  time, 
predicted  a  high-fidelity  readout  of  individual  states  was  possible.  Our  attempts  to  prove  this,  however,  were  limited  by  our  inability,  so 
far,  to  produce  a  pure  excited  state.  Rabi  and  Ramsey  experiments  were  also  both  inconclusive,  because  the  parameters  of  the  boxes 
measured  so  far  were  not  optimal,  requiring  either  very-high  frequency  or  short  duration  pulses,  which  are  difficult  to  realize.  Future 
experiments  in  which  the  fabrication  parameters  of  the  qubit  are  controlled  and  improved  can  greatly  ease  these  requirements,  and  are 
predicted  to  improve  both  the  decoherence  and  relaxation  times  of  the  qubit.  Because  the  backaction  of  our  continuous  SET 
measurement  was  clearly  observable,  we  both  observed  its  dependence  on  various  biasing  parameters  of  the  setup,  and  also 
undertook  theoretical  investigations  of  the  backaction  of  the  SET.  The  initial  steps  on  this  resulted  in  an  invited  “insight  review  article” 
published  in  Nature  in  2000  (see  Appendix).  These  efforts  are  also  continuing  under  the  new  grant  period,  including  an  expanded 
theoretical  collaboration  with  Prof.  Steve  Girvin  at  Yale. 
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Scientific).  Also  appears  on  cond-mat  archive:  cond-mat/0210247. 
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Personnel  and  Budget  Discussion 

Several  postdoctoral  fellows  and  students  were  supported  in  part  or  in  full  by  this  grant,  though  no  student  theses  have  yet 
been  written.  The  main  researcher  supported  by  this  grant  was  a  postdoctoral  fellow,  Dr.  Konrad  Lehnert,  who  performed  many  of  the 
experiments  which  are  reported  here.  Dr.  Lehnert  leaves  the  project  this  December,  to  begin  a  new  position  as  Associate  Fellow  of  the 
Joint  Institute  for  Laboratory  Astrophysics  (JILA)  and  Adjoint  Professor  of  Physics  at  The  University  of  Colorado.  A  second  postdoc,  Dr. 
Ryan  Held,  was  funded  for  a  short  period  of  time  from  this  grant,  concentrating  on  fabrication  and  electromagnetic  simulation  of  qubit 
chips.  The  students  funded  include  Lafe  Spietz,  Ben  Turek,  and  David  Schuster.  Mr.  Turek  and  Mr.  Schuster  now  have  their  tuition  and 
stipend  provided  by  supplementary  SUSPENSE  grants  for  graduate  support,  but  their  equipment  and  supplies  continue  to  be  drawn 
from  this  grant.  A  subcontract  to  Chalmers  University  in  Goteborg,  Sweden  allowed  our  collaborators  there  to  fabricate  the  initial 
devices  with  SETs  and  Cooper-pair  boxes  that  were  used  in  these  experiments.  A  small  amount  of  salary  was  provided  for  Prof  Michel 
Devoret,  during  a  year-long  sabattical  visit  to  our  lab  in  1999-2000.  This  lead  to  a  major  publication  in  Nature,  as  well  as  Prof  Devoret’s 
eventual  move  to  the  US  to  join  the  Yale  Applied  Physics  department.  Some  funds  were  also  used  for  summer  salary  for  the  PI. 

Especially  during  the  one  to  two  years  of  this  project,  there  was  a  major  effort  in  building  apparatuses  and  infrastructure,  and 
therefore  significant  expenses  in  supplies  and  equipment.  A  dilution  refrigerator  was  purchased  with  Pi’s  university  setup  funds,  but 
supplies  for  the  wiring  and  filtering  of  the  refrigerator  for  the  qubit  experiments  were  supplied  by  this  grant,  as  well  as  some  specialty 
equipment  such  as  cryogenic  RF  amplifiers.  Because  of  the  overlap  with  the  renewal  grant  covering  this  project,  funds  from  the  last 
year  of  this  grant  were  rebudgeted  and  applied  towards  the  purchase  of  a  second  dilution  refrigerator,  which  will  double  our 
experimental  capability  when  it  becomes  operation  in  mid-2003.  Other  major  equipment  purchases  included  a  40  GHz  network 
analyzer,  a  superconducting  magnet  and  associated  low-loss  dewar,  microwave  components  for  cabling  and  pulse  generation,  the  fast 
data-acquisition  cards,  finite-element  analysis  software  packages  for  microwave  and  electrostatic  circuit  simulation,  and  RF 
synthesizers  with  fast  pulse  capability  for  generating  the  carrier  for  the  RF-SET  readouts. 
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Detailed  Discussion  of  Technical  Results 


Following  the  general  approach  described  in  the  summary  above,  during  this  funding  period  we  developed  the  capability  for 
performing  single-charge  measurements,  using  the  high-speed  version  of  the  single  electron  transistor  (RF-SET)  as  a  high- 
performance  electrometer.  During  this  project,  we  performed  the  firs  ever  charge  measurements  with  an  RF-SET,  and  have  improved 
the  techniques,  control  software,  and  calibration  procedures  so  that  both  very  high  speed  (submicrosecond)  and  truly  high-precision 
(level  of  a  millielectron  of  charge)  charge  measurements  are  now  routine.  Several  manuscripts  employing  these  results  to  study  the 
single  Cooper-pair  and  single  electron  box  are  now  in  press  or  in  preparation,  and  this  work  will  continue. 

Several  generations  of  chip  layouts  were  fabricated  and  modeled,  and  two  different  designs  have  been  measured  at 
cryogenic  temperatures.  An  atomic  force  microscope  image  of  the  type  of  chip  used  for  the  coherence  measurements,  showing  the 
Cooper-pair  box  (with  a  small  SQUID  loop  to  allow  in-situ  tuning  of  the  Josephson  energy)  and  the  integrated  RF-SET  readout 
electrometer,  is  presented  in  Figure  2.  Two  samples  of  slightly  different  designs  were  studied,  for  which  the  electrometers  had  a 
charge  noise  of  approximately  4  x  10'5  electrons/rt(Hz)  at  high  frequencies. 


Figure  2.  Atomic  force  microscope  image  (left)  of  single  Cooper-pair  qubit  (bottom)  and  RF-SET  measuring  electrometer 
(top),  fabricated  in  AI/AIOx/AI  tunnel  junctions  using  direct-write  electron-beam  lithography.  The  CPB  has  two  junctions  formed  into  a 
superconducting  loop,  which  allows  a  tuning  of  the  effective  Josephson  energy,  Ej,  with  a  magnetic  field.  Right  image  is  results  of  a  3- 
d  finite  element  simulation  of  the  electromagnetics  of  this  structure,  used  to  find  the  full  capacitance  matrix  for  the  circuit. 

In  our  first  measurements  on  the  CPB  qubit,  a  simple  scheme  was  used  in  which  the  electrometer  continuously  measures  the 
time-averaged  charge  on  the  box,  and  acquisition  times  per  point  are  longer  than  the  dynamical  timescales  in  the  box.  This  is  therefore 
an  ensemble  measurement,  with  many  realizations  of  the  single  psuedo-spin  acquired  sequentially  in  time.  Nonetheless,  this  mode 
allows  a  rapid  characterization  of  the  CPB  spectrum.  By  studying  the  average  charge  as  a  function  of  temperature,  several  important 
parameters  of  the  system,  including  the  gate  capacitance,  the  measurement  capacitance  that  couples  the  box  to  the  electrometer,  and 
the  total  box  capacitance  or  charging  energy  Ec,  were  determined. 


Figure  3:  CW  spectroscopy  of  Cooper-pair  box,  as  observed  in  a  continuous  measurement  with  the  RF-SET.  A)  A  high- 
frequency  (35  GHz)  microwave  signal,  plus  a  do  gate  voltage,  are  applied  to  the  CPB.  B)  Diagram  of  energy  levels  vs.  dc  gate  voltage. 
At  two  points,  on  either  side  of  the  degeneracy  point,  the  energy  level  splitting  is  resonant  with  the  microwave  drive,  and  transitions 
between  ground  and  excited  states  are  driven.  On  the  left  hand  side,  the  drive  causes  transitions  from  the  charge  0  to  charge  1  state, 
while  the  opposite  is  true  on  the  right-hand  side,  thus  an  absorption  peak  and  dip,  respectively.  C)  Average  charge  on  the  CPB  as  a 
function  of  normalized  gate  voltage.  Lower  curve  is  without  microwaves,  the  dotted  and  dashed  curves  show  the  expected  ground  and 
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excited  state  charge  signal,  solid  upper  curve  shows  the  resonant  absorption  features  caused  by  the  coupling  of  ground  and  excited 
states.  The  actual  peaks  observed  correspond  to  two-photon  transitions,  where  the  box’s  energy  splitting  is  ~  70  GHz. 


With  this  setup,  the  signal-to-noise  in  the  determination  of  the  charge  state  of  the  box  is  very  high.  In  fact,  when  performing 
CW  spectroscopy  of  the  charge  states  in  the  CPB,  the  photon  resonance  peaks  can  even  be  seen  in  realtime  on  an  oscilloscope.  In 
this  measurement,  a  CW  microwave  excitation  is  applied  to  the  gate  of  the  box,  and  resonant  absorption  peaks  (see  Fig.  3)  are 
observed  where  the  energy  splitting  between  the  ground  and  excited  states  matches  the  applied  photon  frequency.  By  studying  the 
photon  peak  height  and  width  as  a  function  of  the  microwave  power  applied,  we  can  extract  (see  Fig.  4)  the  ensemble  averaged 
decoherence  time,  T2*,  using  standard  results  for  the  driven  two-level  system.  The  measured  coherence  time  of  about  1  nanosecond  is 
similar  to  that  observed  without  echo  methods  by  Nakamura,  and  probably  limited  by  the  low-frequency  1/f  offset  charge  noise  of  the 
device.  However,  our  initial  use  of  a  continuous  measurement  means  that  additional  dephasing  due  to  the  SET  electrometer  may  be 
present.  A  determination  of  the  dominant  dephasing  mechanism,  and  a  separation  of  the  measurement  effects,  are  now  underway  and 
will  be  extended  in  future  work.  We  have  even  observed  that  the  SET  is  able  to  drive  transitions  (“state  mixing”)  between  states  in  the 
box,  which  limits  the  investigations  to  large  energy  splittings,  away  from  the  degeneracy  point  in  gate  voltage.  Using  time-gated 
measurements  with  the  RF-SET,  we  have  measured  this  mixing  time  due  to  the  SET  is  about  10  microseconds. 
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Figure  4:  A)  Expanded  view  of  a  photon  peak  from  Figure  3.  B)  The  power  dependence  of  the  linewidth  of  the  photon  peaks 
in  Figure  3.  The  square  of  the  observed  photon  peak  width  is  plotted  against  the  microwave  power,  or  Rabi  frequency  squared.  The 
linear  behavior  allows  the  extraction  of  a  T2*  value  for  the  CPB  of  ~  1  nanosecond. 

One  of  the  first  tasks  that  we  needed  to  perform  with  this  system  is  to  determine  the  parameters  of  our  box  and  electrometer 
circuit,  and  to  verify  that  the  Hamiltonian  of  our  circuit  matches  that  of  a  driven  two-level  or  qubit.  The  expected  Hamiltonian  of  the 
Cooper-pair  box  circuit  is 


H  =  —2EC  (1  -  2ng  )oz 


(Equation. 1), 


where  <TX  andoz  are  Pauli  spin  matrices  and  ng  is  total  polarization  charge  applied  to  the  gate  electrode, 
n:,  =CgVg/2e-noff  ,  in  units  of  a  Cooper-pair's  charge.  The  offset  charge  noff  accounts  for  the  uncontrolled  potential  from  charges 
nearby  the  box  island.  The  charging  energy,  Ec  =e2 HCZ ,  is  the  electrostatic  energy  required  to  add  one  electron  to  the  island  and  the 
Josephson  energy,  £'™x  =  Ti nEI\e1RJ ,  is  the  effective  tunneling  matrix  element  for  Cooper-pairs  across  a  junction  with  resistance  R, 
in  a  superconductor  with  BCS  gap  A  .  The  junction  is  in  fact  a  composite  of  two  junctions  shorted  together  to  form  a  loop  with 
1  (pm)2  area  (Fig.  2).  The  effective  Josephson  energy  E ;  of  the  pair  of  junctions  is  tuned  by  introducing  magnetic  flux  fl>  into  this  loop, 
as  Ej  =  cosfrc®  / ®0) ,  where  ®0  is  the  quantum  of  flux  (h/2e) . 

This  Hamiltonian  was  verified  by  performing  the  CW  spectroscopy  with  continuous  measurement  (as  described  above),  by 
applying  a  continuous  microwave  stimulus  to  the  gate  of  the  Cooper-pair  box,  and  tuning  the  parameters  of  the  TLS  to  find  the 
resonance  condition.  In  Fig.  3B  we  show  the  ground-state  and  excited-state  energies  of  the  Hamiltonian  of  Eq.  1  as  a  function  of  the 
reduced  gate  voltage,  ng .  Figure  3B  shows  that  as  n  is  tuned  through  1/2,  the  charge  of  the  ground  and  excited  states  exchange, 


with  the  ground  state  acquiring  a  charge  of  one  Cooper  pair  and  the  excited  state  attaining  zero  charge.  Also  plotted  is  an  ensemble- 
averaged  measurement  of  the  charge  on  the  box  as  a  function  of  ng  .  It  is  apparent  that  the  box  does  not  remain  in  its  ground  state 

over  a  range  in  0.3<  ng  <0.7.  This  behavior  is  caused  by  backaction  generated  by  currents  flowing  through  the  RF-SET.  We  proceed  by 
studying  the  box  in  the  range  of  ng  where  it  does  remain  in  its  ground  state.  When  a  38  GHz  microwave  signal  is  applied  to  the  gate, 
we  observe  clear  evidence  that  the  box  is  a  coherent  two-level  system.  Resonant  peaks  appear  (Fig.  3C)  in  Qbox  that  are  sharp  and 
symmetrically  spaced  about  n  =  0.5  .  The  two  features,  a  peak  for  ng  <  0.5  and  a  dip  for  ng  >  0.5  ,  both  correspond  to  the  change  in 
Qbox  when  the  box  spends  some  time  in  the  excited  state.  Because  Qhax  is  an  average  of  thousands  of  repeated  measurements,  the 
peak  height  indicates  the  probability  of  finding  the  box  in  its  excited  state. 

The  resonant  peaks  permit  a  spectroscopic  determination  of  the  parameters  in  the  box's  Hamiltonian.  By  plotting  the  locations 
of  the  resonances  as  a  function  of  ng  and  <3> ,  we  find  good  agreement  between  this  spectroscopy  and  the  difference  Em  =  hcom 

between  ground-state  and  excited-state  energy  expected  from  Eq.  1.  An  independent  measurement  of  Ec  ,  observing  the  thermal 
smearing  of  the  Coulomb  staircase  in  the  normal  state,  demonstrates  that  these  peaks  occur  when  the  irradiating  frequency  co  is  half 
com  ,  indicating  that  these  peaks  correspond  to  a  2-photon  transition.  We  find  a  single  value  for  Ec  and  for  Zs“x  that  account  for  the 
location  of  the  resonant  peaks  at  applied  frequencies  between  32  and  38  GHz  giving  resonant  peaks  for  com  between  64  and  76  GHz 
(Fig.  5A).  At  lower  frequencies,  the  2-photon  peaks  would  appear  at  an  ng  for  which  the  box  does  not  stay  in  the  ground  state  while 
being  measured  and  are  therefore  not  visible  (For  lower  frequencies,  we  have  observed  the  3-photon  peaks  whose  location  in  n  is 
consistent  with  our  interpretation).  We  are  able  to  extract  the  parameters  of  the  Hamiltonian,  4Ec/h  =149.1  ±0.4  GHz  and 
E™*  lh  =  13.04  +  0.2  GHz  .  The  uncertainties  arise  from  the  systematic  deviation  from  linearity  of  the  function  generator,  which  created 
the  ramp  voltage  for  Vg .  Because  these  measurements  were  made  at  a  temperature  T  =  0.015  K  ,  they  are  in  the  limit 
kj  «  EJ"  «  4  Ec . 


Consistent  with  the  behavior  of  a  TLS,  the  peaks  disappear  for  <J>  =  <I>0/2when  Ej  approaches  0,  which  demonstrates 
that£j  provides  the  coupling  between  the  charge  states  (Fig.  5B).  This  can  be  understood  geometrically  from  the  fact  that  an 
oscillating  gate  voltage  with  amplitude  V“e  =2 en“/C  adds  a  term  to  the  Hamiltonian  in  Eq.  1  nj  cos((ot)aI ,  which  is  collinear  with 
ground  state  of  the  quasi-spin  described  by  Eq.  1  when  Ej  =  0 .  The  microwave  excitation  therefore  applies  no  torque  which  could 
excite  the  quasi-spin  from  its  ground  state. 


(Wo) 


Figure  5  A.)  The  symbols  (open  circles)  represent  the  locations  of  resonant  peaks  as  a  function  the  two  control  parameters 
of  the  Hamiltonian,  ng  and  <5  ,  for  three  frequencies,  co  =  32  (top),  35,  and  38  (bottom)  GHz.  The  three  lines  are  a  fit  to  the  data,  using 

Eq.  1  to  find  n,,  given  com  and  <J>  for  (om  =  64,  70  and  76  GHz.  B.)  The  size,  in  electrons,  of  a  38  GHz  resonant  peak  as  a  function  of 
<f>  (squares).  Note  that  the  peaks  disappear  when  Ej  approaches  zero. 
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In  order  to  measure  the  excited-state  lifetime  7) ,  we  excite  the  box  and  then  measure  the  time  required  to  relax  back  to  the 
ground  state.  A  38  GHz  signal  is  continuously  applied  to  the  gate  and  the  box  gate  is  tuned  to  n  =  0.248  and  E}  =  E™x  so  that  the 
microwaves  resonantly  couple  the  ground  and  excited  state  through  a  2-photon  transition.  Abruptly,  n,,  is  then  shifted  to  n,,  =0.171  in 
30  ns,  slowly  enough  to  be  adiabatic  but  much  faster  than  7) .  The  microwave  excitation  no  longer  resonantly  couples  the  ground  and 
excited  state,  and  the  probability  of  being  in  the  excited  state  decays  in  a  time  7, .  By  averaging  many  of  the  transient  responses  to  this 
stimulus,  we  find  Tt  =1.3  jus  (Fig.  6).  The  lifetime  is  a  quantity  which  is  insensitive  to  slow  drifts  in  noff  and  demonstrates  that  the 
intrinsic  quality  factor  of  the  TLS,  0,  =  Ttcom  =5xl05 . 


time  (as) 


Figure  6:  The  average  charge  on  the  box  is  plotted  (triangles)  as  a  function  of  time  f,  relative  to  t  =  0,  when  ng  is  shifted  from 
ng  -  0.248  to  0.171  in  30  ns.  With  38  GHz  microwaves  applied,  the  shift  in  gate  brings  the  box  out  of  resonance  with  the  microwave 
excitation.  The  decay  in  the  charge  on  the  box  represents  the  decay  of  the  probability  of  the  box  being  in  its  excited  state  after  the  shift 
in  n  .  An  exponential  fit  to  the  data  implies  Tx  =1.3  jus  (line). 


We  can  compare  this  long  lifetime,  with  the  spontaneous  emission  rate  induced  by  the  quantum  fluctuations  of  a  generic 
electromagnetic  environment.  Calculating  the  spontaneous  emission  rate  using  Fermi's  golden  rule  gives 


J_ 
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^ g 
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-  sin  (8)Sr(co01) 
n , 


(2) 


where  Sy(co)  =  27?<w(Re(Z0))  is  the  voltage  spectral  density  (per  Hz)  of  the  quantum  fluctuations  of  an  environment  with  an  impedance 

Z0  at  frequency  co  and  sin  9  =  Ej  /hco0x .  The  quantity  CTg  is  the  total  capacitance  of  the  box  to  nearby  metal  traces,  including 

intentional  coupling  to  the  gate  lead  and  other  unintended  capacitive  coupling  (Fig.  1).  We  calculate  7j  fora  Z0=50Q  environment  to 

be  between  1  and  4  jus  ,  extracting  CTg  =  60  aF  with  about  50%  uncertainty  from  an  electrostatic  simulation  of  the  chip  layout.  We  do 

not  claim  to  have  demonstrated  that  the  lifetime  is  limited  by  spontaneous  emission;  however,  if  there  are  additional  relaxation 
processes,  either  due  to  the  nearby  electrometer  or  fluctuations  of  some  microscopic  degree  of  freedom  in  the  box,  their  influence  is  at 
most  comparable  to  that  of  spontaneous  emission  into  a  50  Q.  environment.  This  very  long  lifetime  emphasizes  the  importance  of 
understanding  and  reducing  the  excess  dephasing. 


Our  measurement  of  the  excited  state  lifetime  of  box  is  remarkable  for  two  reasons.  First,  it  shows  that  a  quantum-coherent 
microelectronic  circuit  can  have  a  lifetime  that  approaches  the  limit  set  by  spontaneous  emission  of  a  photon  into  the  electromagnetic 
environment.  Furthermore,  it  is  achieved  while  the  two-level  system  is  continuously  measured.  This  means  that  the  coherence  time  in 
the  Cooper-pair  box  can  be  long  lived,  if  the  sources  of  inhomogeneous  decoherence  can  be  reduced.  M.  Devoret  and  coworkers 
(formerly  at  Saclay,  and  now  at  Yale)  have  indeed  observed  dephasing  times  close  to  the  inelastic  lifetime  limit  in  a  similar  system,  by 
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operating  at  the  charge  degeneracy  point,  ng  =  0.5  .  Furthermore,  given  the  observed  electrometer  sensitivity  of  4x10  5  e/VHz  ,  the 

excited-state  lifetime  is  long  enough  that  a  single  measurement  can  discriminate  between  the  box  in  its  excited  state  and  the  box  in  its 
ground  state  with  99%  confidence.  Both  of  these  are  vital  to  implementing  a  quantum  computer.  In  future  work  we  intend  to  operate 
these  CPB  qubits  at  a  charge-noise  insensitive  point  in  gate  bias  (ng=0.5),  and  also  to  perform  single-shot  readouts  of  the  quantum 
state,  both  of  which  are  required  before  proceeding  to  studies  of  entanglement  in  superconducting  qubits  and  implementation  of 
quantum  logic  gates. 

One  attractive  feature  of  the  RF-SET  measurements  we  perform  is  that  the  operating  point  of  the  SET,  and  its  backaction 
effects  on  the  CPB,  can  be  widely  varied.  In  fact,  we  have  already  observed  a  dependence  of  the  average  charge  state  of  the  CPB,  i.e. 
the  degree  of  polarization  of  our  psuedo-spin,  on  this  operating  point  (see  Fig  9).  We  find  that  higher  bias,  and  higher  drain-source 
current,  lead  to  an  increase  in  “mixing”,  or  driven  transitions  between  states  in  the  CPB.  We  could  further  observe  the  evolution  of  the 
quantum  state  of  the  CPB  as  a  function  of  time  after  turning  on  the  measurement,  which  allowed  us  to  determine  a  “mixing  time”  for 
the  states  of  about  10  ps.  By  repeating  and  refining  these  types  of  measurements,  and  by  tuning  the  Josephson  energy  and  energy 
level  splitting  of  the  box,  we  can  use  the  CPB  two-level  system  as  a  spectrum  analyzer  to  measure  the  full  noise  spectrum  produced  by 
the  SET.  Another  interesting  possibility  is  to  observe  the  dephasing  effects  of  the  measurement  on  the  width  of  photon  resonances,  i.e. 
the  actual  dephasing. 
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Figure  7:  Dependence  of  CPB  states  on  bias  conditions  of  the  SET  electrometer.  The  top  panel  shows  the  conductance  through 
the  measuring  SET  as  a  function  of  the  voltage  on  the  SET  drain  source  (VdS,  vertical  axis)  and  gate  (Vge,  horizontal  axis)  voltage.  The 
arrows  indicate  several  operating  points  which  have  been  investigated.  Bottom  panel  shows  the  average  charge  on  the  CPB  as  a 
function  of  box  gate  voltage  that  is  observed  for  each  operating  point. 

The  RF-SET  is  an  appealing  measurement  device  not  only  because  of  its  tunability,  but  because  the  measurement  process  itself 
can  be  understood  quantum  mechanically.  Most  theoretical  treatments  of  this  problem,  however,  have  focused  on  the  simpler  case  of 
a  normal  SET.  In  this  limit,  it  is  expected  that  the  SET  can  be  nearly  quantum  limited,  i.e.  the  dephasing  introduced  by  the 
measurement  approaches  that  required  by  the  uncertainty  principle.  However,  the  actual  measuring  SET  is  superconducting,  leading 
to  more  complicated  dynamics  and  operating  possibilities  (e.g.  the  complicated  charging  diamonds  of  Fig.  7).  In  our  first  samples,  we 
found  that  the  least  invasive  continuous  measurement  was  obtained  at  low  drain-source  bias  on  the  SET,  using  a  resonant  Cooper- 
pair/quasiparticle  tunneling  feature  (or  “double  JQP”).  We  are  further  exploring  these  types  of  features  in  the  superconducting  SET, 
and  have  begun  to  compare  with  theoretical  models  developed  in  collaboration  with  Steve  Girvin’s  group  at  Yale.  We  will  explore  the 
various  parameters  of  the  SET  such  as  resistance,  capacitance,  and  coupling  to  the  CPB,  and  how  to  optimize  these  for  both 
continuous  and  pulsed  quantum  measurements  in  future. 


Report  of  Inventions 

No  inventions,  patentable  or  otherwise,  to  report. 

Appendix:  Copies  of  Major  Publications 

Three  publications,  as  described  above,  are  appended  to  this  electronic  file. 
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Transistors  have  continuously  reduced  in  size  and  increased  in  switching  speed  since  their  invention  in  1947. 
The  exponential  pace  of  transistor  evolution  has  led  to  a  revolution  in  information  acquisition,  processing  and 
communication  technologies.  And  reigning  over  most  digital  applications  is  a  single  device  structure  —  the 
field-effect  transistor  (FET).  But  as  device  dimensions  approach  the  nanometre  scale,  quantum  effects 
become  increasingly  important  for  device  operation,  and  conceptually  new  transistor  structures  may  need  to 
be  adopted.  A  notable  example  of  such  a  structure  is  the  single-electron  transistor,  or  SET1-4.  Although  it  is 
unlikely  that  SETs  will  replace  FETs  in  conventional  electronics,  they  should  prove  useful  in  ultra-low-noise 
analog  applications.  Moreover,  because  it  is  not  affected  by  the  same  technological  limitations  as  the  FET, 
the  SET  can  approach  closely  the  quantum  limit  of  sensitivity.  It  might  also  be  a  useful  read-out  device  for  a 
solid-state  quantum  computer. 


It  is  now  possible  to  put  a  billion  transistors  on  a 
single  chip  operating  with  a  clock  period  of  a 
billionth  of  a  second.  Most  probably,  the  trend  in 
reducing  dimensions  and  times  will  continue  in  the 
next  decade.  But  as  the  number  and  density  of  gates 
and  memory  elements  increase,  the  energy  of  signals  also 
has  to  be  reduced  to  keep  the  power  dissipation 
sufficiently  low. 

Surprisingly,  even  though  the  size  of  a  typical  transistor 
in  a  microcomputer  chip  is  now  j  ust  a  few  hundred  nanome¬ 
ters,  its  functioning  remains  essentially  classical:  quantum 
mechanics  only  enters  in  the  explanation  of  the  values  of  the 
physical  parameters  of  materials,  like  the  band-gap  of  a 
semiconductor.  Otherwise,  the  discreetness  of  matter  and 
the  wave-like  properties  of  electrons  can  be  largely  ignored 
in  the  understanding  of  the  behaviour  of  electrical  signals  in 
today’s  integrated  circuits. 

But  as  devices  get  smaller,  faster  and  more  densely  packed, 
quantum  effects  will  have  increasingly  to  be  taken  into 
account.  Even  well  before  we  reach  the  ultimate  limit  where 
transistors  are  reduced  to  the  size  of  an  atom  or  a  molecule, 
we  encounter  four  limits.  Quantum  phenomena  become 
significant  when  (1)  signal  energy,  (2)  signal  charge,  (3) 
device  dimension,  and  (4)  device  size  tolerance  approach, 
respectively,  the  energy  of  one  photon,  the  charge  of  one 
electron,  the  electron  wavelength,  and  the  size  of  one  atom. 

Much  research  has  been  devoted  to  assess  if  quantum 
effects  arising  from  these  conditions  will  force  the  adoption 
of  new  physical  principles  or  if  they  can  simply  be  tamed  by 
better  control  of  the  chip  structure  at  the  atomic  level.  To  our 
knowledge,  there  is  no  general  consensus  on  the  answer  to 
this  question. 

Another  research  direction  has  been  to  exploit  quantum 
effects  arising  in  devices  of  nanometer  scale  to  implement  a 
function  that  cannot  be  performed  by  present  devices.  In 
some  applications,  which  operate  at  limits  ( 1 )  or  ( 2)  or  both, 
it  is  not  only  inevitable  but  also  desirable.  In  astronomy,  for 
instance,  it  is  important  to  extract  as  much  information  as 
possible  from  a  single  photon5.  Recent  advances  in  quantum 
information  theory6  indicate  that  scalable  switching 
elements  that  behave  fully  as  quantum  systems  would  not 


simply  make  calculations  with  minimal  energy,  they  could 
in  addition  perform  tasks  that  would  be  impossible  with 
conventional  computers.  In  a  quantum  computer,  usual  bits 
are  replaced  by  quantum  bits  or  ‘qubits’  which  can  be 
‘entangled’  with  each  other,  thereby  carrying  a  new  type  of 
information  that  is  useful  in  solving  highly  parallel  tasks. 
New  types  of  devices  are  needed  to  read-out  such  qubits, 
that  is,  to  amplify  their  associated  single-quantum  signals. 

In  the  realm  of  atomic  physics  and  quantum  optics,  the 
detection  of  individual  microscopic  particles  travelling  in 
vacuum,  such  as  photons  and  electrons,  is  now  performed 
routinely  with  almost  unity  efficiency  by  instruments 
derived  from  the  photomultiplier.  However,  the  measure¬ 
ment  of  electrical  signals  resulting  from  the  motion  of  a 
single  electron  in  a  circuit  involves  an  amplifier  having  not 
only  a  good  energy  sensitivity,  but  also  electrical  characteris¬ 
tics  that  are  adapted  to  this  circuit. 

A  particularly  simple  and  spectacular  example  of  such  a 
device  is  the  single-electron  transistor  (SET)1-3,  which  exploits 
the  quantum  phenomenon  of  tunnelling  to  control  and 
measure  the  movement  of  single  electrons  inside  a  solid-state 
circuit.  SETs  are  extremely  precise  solid-state  electrometers4'7, 
already  out-performing  state-of-the-art  conventional  transis¬ 
tors8  by  three  orders  of  magnitude.  Their  charge  sensitivity  has 
been  shown  to  be  as  low  as  a  few  10-5  e/VHz,  which  means  that 
a  charge  variation  of  1 0  -  5  e  can  be  detected  in  a  measurement 
time  of  1  s  (the  precision  improves  as  the  square  root  of  the 
measurement  time).  This  result  is  only  an  order  of  magnitude 
away  from  the  theoretical  limit  of  10-6  e  /  VHz.  SETs  have 
applications  in  metrology9  and  single-photon  detection10,11. 
Furthermore,  it  has  been  realized  in  the  past  fewyears  that  SETs 
can  perform  a  measurement  on  a  single  quantum  two-state 
system  (qubit) 1213,  perturbing  its  quantum  evolution  in  a  mini¬ 
mal  way.  That  is,  the  SET  is  a  charge  amplifier  operating  in  the 
vicinity  of  the  quantum  limit.  It  would  be  a  practical  read-out 
device  for  several  solid-state  implementations  of  qubits14.  In 
this  article  we  review  these  latest  developments. 

Conventional  transistors 

Before  discussing  the  performance  of  SETs,  it  is  useful  to  recall 
the  operating  principle  of  the  most  common  transistor, 
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the  metal-oxide  semiconductor  field-effect  transistor  (MOSFET). 
Figure  1  depicts  schematically  the  layout  of  the  device  and  its  operating 
principle  (we  restrict  ourselves  here  to  the  nMOSFET,  in  which  the 
majority  charge  carriers  are  electrons).  Two  conducting  electrodes, 
called  the  source  and  drain,  are  connected  by  a  channel  made  of  a 
material  in  which  the  number  of  conduction  electrons  can  be  varied,  in 
practice  a  semiconductor  (Fig.  la).  A  voltage  is  applied  to  the  ‘gate’,  a 
third  conducting  electrode  that  is  separated  from  the  channel  by  a  thin 
insulating  layer.  When  the  gate-source  voltage  is  zero,  there  are  no 
conduction  electrons  in  the  channel  as  the  effective  potential  they 
would  experience  there  is  larger  than  in  the  leads  (Fig.  lb).  The  channel 
is  therefore  in  an  insulating  state.  But  when  the  gate  voltage  is  increased 
with  respect  to  the  source,  the  potential  experienced  by  the  electrons  in 
the  channel  decreases  and  they  populate  the  channel  just  under  the  gate 
(Fig.  1  c) .  The  channel  becomes  conducting.  The  larger  the  gate  voltage, 
the  larger  will  be  the  channel  electron  population  that  can  participate 
in  the  current.  Eventually,  all  the  electron  states  in  the  energy  window 
set  by  the  source-drain  voltage  can  propagate  through  the  channel  and 
the  current  no  longer  depends  on  the  gate  voltage.  This  saturation 
regime  is  depicted  in  Fig.  Id. 

This  field  effect  provides  an  amplification  mechanism  since  an 
increase  in  gate  voltage,  bringing  a  modest  current  to  the  gate 
electrode,  can  switch  on  a  larger  current  through  the  channel  (Fig.  2). 
The  source-drain  current  is  determined  by  the  conductance  of  the 
channel,  which  in  turn  depends  on  two  factors:  the  density  of  its 
conduction  electrons  and  their  mobility.  The  electron  density  is 
controlled  directly  by  the  gate  voltage.  The  electron  mobility  is  set  by 
the  collisions  of  electrons  with  static  irregularities  of  the  crystal  as 
well  as  with  its  dynamic  deformations  due  to  thermal  agitation. 
When  thermal  agitation  is  the  predominant  factor,  electron  mobility 
is,  to  a  large  extent,  independent  of  the  gate  voltage.  However,  at  low 
temperatures,  mobility  can  also  increase  when  the  density  increases, 
reinforcing  the  influence  of  the  gate  voltage  on  the  source-drain 
current. 

Note  that  so  far  we  have  made  no  reference  to  the  wave-like 
properties  of  electrons  nor  to  the  fact  that  the  channel  is  made  from 
individual  atoms.  The  only  quantum  property  that  has  had  a  role  in 
our  explanation  is  the  Pauli  principle,  which  dictates  that  each  possi¬ 
ble  state  for  an  electron  in  the  channel  can  be  occupied  at  most  by 
only  one  electron.  This  means  that  only  a  certain  number  of  electrons 
can  accumulate  in  the  channel,  setting  a  limit  on  the  current  flow. 

However,  the  quantum  properties  of  electrons  and  atoms  will  be 
increasingly  important  as  FETs  are  made  smaller.  For  example,  the 
wave  nature  of  electrons  will  influence  the  way  they  travel  through 
the  channel.  When  the  transverse  dimension  of  the  channel  becomes 
comparable  to  the  wavelength  of  electrons  (around  100  nm), 
electron  propagation  becomes  more  sensitive  to  the  atomic  disorder 
in  the  device,  which  is  inherent  in  the  present  fabrication  process. 
The  disorder  makes  the  channel  remain  insulating  even  when  the 
density  of  electrons  is  increased.  Effects  of  this  kind,  which  result 
from  reaching  the  limit  to  device  dimension  (limit  (3)  above),  pose  a 
major  problem  if  the  reduction  of  size  is  not  accompanied  by  an 
improvement  in  the  atomic  structure  of  the  fabricated  devices. 

If,  however,  the  atomic  structure  of  the  FET  could  be  made 
defect-free,  a  recent  analysis15  shows  it  would  continue  to  function  in 
the  regime  where  the  electron  propagation  in  the  channel  is  wave¬ 
like.  It  would  thus  ‘break’  limit  (3)  and  could  be  further  scaled  down 
until  a  channel  length  of  about  8  nm  is  reached,  a  stage  at  which  limit 
(4)  seems  to  severely  affect  the  performance  of  the  FET. 

In  an  ideal  situation,  we  might  want  electrons  to  be  scattered  only 
by  the  gate-dependent  potential,  with  other  scattering  mechanisms 
always  being  detrimental  to  amplification.  But  the  confinement  of 
electrons  in  the  channel  by  tunnel  barriers,  in  the  source-drain 
direction,  can  also  lead  to  a  new  kind  of  amplification  principle.  This 
is  the  basis  of  the  SET.  As  we  shall  see,  this  new  principle  circumvents 
the  problem  of  the  FET  that  the  gate  capacitance,  which  is  an  impor¬ 
tant  factor  in  determining  charge  sensitivity,  is  tied  to  the  size  of  the 
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Figure  1  Principle  of  a  metal-oxide  semiconductor  field-effect  transistor  (MOSFET). 
a,  The  device  consists  of  two  conducting  electrodes  (source  and  drain)  connected  by  a 
semiconducting  channel.  The  channel  is  influenced  electrostatically  by  the  gate,  a 
conducting  electrode  separated  from  the  channel  by  thin  insulating  layer,  b,  When  the 
voltage  applied  between  the  gate  and  the  source  is  zero,  the  Fermi  energy  of  the  source 
and  drain  lies  in  the  gap  of  the  semiconductor.  Here,  we  sketch  the  potential  (red  curve) 
seen  by  conduction  electrons  when  they  travel  from  source  to  drain  along  a  line  in  the 
channel  just  under  the  gate.  There  are  no  filled  electron  states  (green  lines)  in  the 
channel,  which  as  a  result  remains  insulating,  c,  When  the  gate  voltage  is  increased, 
the  potential  seen  by  conduction  electrons  is  lowered.  There  are  now  filled  states  in  the 
channel  at  the  Fermi  energy  of  the  source  and  drain.  The  device  conducts,  d,  When  the 
gate  voltage  is  increased  further,  the  current  finally  saturates  when  all  the  states  in  the 
bias  window  are  filled. 


channel.  Until  the  technology  for  reliably  fabricating  FETs  with 
channels  of  nanometre-scale  dimensions  has  been  developed,  the 
SET  is  the  best  device  in  terms  of  charge  sensitivity. 

Operating  principle  of  SETs 

Unlike  the  FET,  whose  principle  does  not  require  the  motion  of  elec¬ 
trons  to  be  quantum-mechanical,  the  SET  is  based  on  an  intrinsically 
quantum  phenomenon:  the  tunnel  effect  through  a  metal- 
insulator-metal  junction.  When  two  metallic  electrodes  are  separated 
by  an  insulating  barrier  whose  thickness  is  only  ~  1  nm,  electrons  at  the 
Fermi  energy  can  traverse  the  insulator  even  though  their  energy  is  too 
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Figure  2  Variation  of  the  source-drain  current  in  a  MOSFET  as  a  function  of  the  gate 
voltage.  When  the  gate  voltage  is  increased  from  zero,  the  source-drain  current  is 
turned  on.  This  device  can  be  used  both  in  digital  electronics  and  as  an  amplifier  for 
analog  signals. 


low  to  overcome,  in  a  classical  motion,  the  large  potential  barrier  of  the 
insulating  region.  The  tunnel  effect  manifests  itself  by  a  finite 
resistance  Rr  of  the  insulating  barrier.  This  resistance  depends  both  on 
the  transmission  coefficient  9~  of  the  barrier  to  electron  waves  (which 
is  an  exponentially  decreasing  function  of  its  thickness)  and  on  the 
number  M  of  independent  electron  wave  modes  impinging  on 
the  barrier  (this  number  is  equal  to  the  area  of  the  junction  divided 
by  the  square  of  the  electron  wavelength) .  The  SET  uses  a  key  property 
of  the  tunnel  effect  in  a  many-electron  system:  for  barriers  such  that 
2TM  1,  the  charge  Q  transferred  through  the  barrier  becomes 
quantized  with  Q  =  Ne,  where  N  is  an  integer16.  In  other  words,  for  N 
not  to  be  subject  to  quantum  fluctuations,  the  resistance  of  the  junc¬ 
tion  must  be  large  compared  with  the  resistance  quantum  RT  hle2  = 
_RK  =  25.8  kfl  (refs  17,18). 

The  SET  consists  of  two  such  tunnel  junctions  placed  in  series 
(Fig.  3a, b).  An  ‘island’  is  thus  formed  between  the  two  junctions.  A 
gate  electrode  is  coupled  electrostatically  to  the  island.  The  SET  can 
thus  be  described  as  a  FET  in  which  the  semiconducting  channel  has 
been  replaced  by  a  metallic  island  sandwiched  between  two  tunnel 
barriers.  The  island  has  a  total  capacitance  CE,  which  the  sum  of  the 
gate  and  junction  capacitances  Cx  =  Cg  +  Cn  +  CJ2. 

If  the  dimensions  of  the  island  are  sufficiently  small,  the  charging 
energy  Ec  =  e2/(2Cx)  of  one  extra  electron  in  the  island  will  become 
larger  than  the  energy  of  thermal  fluctuations:  Ec  >  kB  T,  where  fcB  is 
the  Boltzmann  constant  and  Tis  temperature.  In  practice,  for  devices 
fabricated  by  standard  electron-beam  lithography,  Cx  is  of  the  order 
of  a  femtofarad  and  the  charging  energy  is  of  order  1  K,  necessitating 
temperatures  below  300  mK  to  satisfy  the  above  charging  energy 
criterion.  Over  the  past  few  years,  however,  experiments  have 
shown  that  with  advanced  fabrication  methods,  room  temperature 
operation  is  possible19-21. 

Because  electrons  interact  strongly  via  the  Coulomb  interaction 
when  they  pass  through  the  island,  the  analysis  of  the  SET  differs 
fundamentally  from  that  of  the  FET.  In  the  FET,  electrons  go  from 
source  to  drain  independently,  and  in  such  numbers  that  one  can 
consider  that  the  potential  seen  by  one  is  an  average  which  does  not 
depend  on  the  configuration  adopted  by  all  of  the  others.  Electrical 
transport  results  from  a  simple  addition  of  the  motion  of  each 
electron.  In  the  SET,  by  contrast,  transport  results  from  transitions 
between  collective  charge  states  of  the  system.  These  charge  states  are 
described  by  the  two  numbers  Nr  and  N2  of  electrons  having  traversed 
thejunctions  (Fig.  3b). 

The  behaviour  of  the  device  is  governed  by  the  global  electrostatic 
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Figure  3  The  single-electron  tunnelling  transistor  (SET),  a,  Simplified  three- 
dimensional  structure  of  the  SET.  The  channel  of  the  FET  is  replaced  here  by  a 
sandwich  consisting  of  a  nanoscale  metal  electrode  (island),  which  is  connected  to  the 
drain  and  the  source  by  tunnel  junctions.  As  in  the  FET,  a  gate  electrode  influences  the 
island  electrostatically,  b,  Circuit  diagram  of  the  SET.  The  square  box  symbol  represents 
a  tunnel  junction,  and  integers  /V,  and  N2  denote  the  numbers  of  electrons  having 
tunnelled  through  the  two  junctions.  Each  junction  is  characterized  by  its  capacitance 
and  its  tunnel  resistance. 


energy16  £el  =  EC[N2  -  N,  -  (CgVg/e)  -  (C2VJe)  +  q0]2  -  eN2Vis, 
which  includes  the  energy  stored  in  the  junction  and  gate  capaci¬ 
tances,  as  well  as  the  work  done  by  the  voltage  sources.  Here,  V  and 
Vis  are  the  voltages  applied  between  gate  and  source,  and  drain  and 
source,  respectively.  The  so-called  offset  gate  charge  q0  is  a  phenome¬ 
nological  quantity  describing  the  fact  that  electric  fields  in  the 
capacitances  of  the  system  are  non-zero  even  when  the  island  is 
neutral  and  when  no  voltage  is  applied.  It  takes  a  randomly  different, 
non-integer  value  for  each  device  and  cool-down.  It  also  fluctuates 
slowly  in  time  with  a  1// spectral  density22.  We  will  discuss  its  effect  in 
more  detail  below.  But  as  far  as  the  amplification  mechanism  of  the 
SET  is  concerned,  we  can  treat  it  as  a  constant. 

According  to  the  so-called  ‘global  rules’,  also  known  as  ‘orthodox 
theory’,  tunnel  events  will  take  place  independently  on  each  junction 
at  a  rate  governed  by  the  global  energy,  provided  that  the  junction 
resistances  satisfy  RT1,  RT2  >  RK  and  that  the  voltage  sources  Vg  and  Vds 
have  negligible  internal  impedance,  on  the  scale  of  the  resistance 
quantum,  around  the  Coulomb  frequency  Eclh  (ref.  17). 

In  this  regime,  each  tunnel  event  creates  one  electron-hole  pair, 
the  electron  and  the  hole  being  on  opposite  sides  of  the  j  unction.  The 
succession  of  tunnel  events  constitutes  a  Poisson  process.  More 
specifically,  a  tunnel  event  will  take  place  on  junction  i  with  a  rate 
given  by  =  [H(RTie2)}  [A£j/(1  -  exp(  -  A£,/fcBT))]  where  A£,-  = 
£el{N“ore,  N, }  -  £el{N,after,  Nj}. 

At  zero  temperature,  tunnel  events  take  place  only  if  they  are  ener¬ 
getically  allowed,  that  is,  A£,  >  0.  For  a  drain-source  voltage  below  the 
Coulomb  gap  voltage  e/Cx,  the  current  therefore  depends  critically 
on  the  value  of  the  gate  voltage.  If  the  gate  voltage  is  such  that 
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a  Cg  l/g  -  eqg  =  0  mod  e 
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b  Cg  l/g  -  eqg  =  e/2  mod  e 


p  p+1  p+2  p+3  (A/-|+A/2)/2 

Figure  4  Charge  levels  of  the  SET.  a,  Charge  levels  in  the  full  Coulomb  blockade; 
b,  charge  levels  in  the  lifted  Coulomb  blockade  case.  The  integers  A/,  and  l\l2  denote  the 
numbers  of  electrons  having  tunnelled  through  the  two  junctions  (see  Fig.  3).  In  a,  the 
only  possible  conduction  process  is  co-tunnelling,  that  is,  a  weak  second-order  process 
involving  the  tunnelling  of  two  electrons  on  two  different  junctions  and  a  virtual  charge 
state  (dashed  purple  arrows).  In  b,  the  main  conduction  process  is  a  cascade  of  tunnel 
events  involving  the  two  junctions  sequentially  (full  purple  arrows). 


CgVg  =  q0e  mod  e,  the  island  has  a  well  defined  number  of  electrons. 
Tunnel  events  cannot  take  place  on  the  junctions  because  the  global 
energy  would  increase  (Fig.  4a)  and  the  current  is  zero.  This  is  the 
phenomenon  of  Coulomb  blockade.  But  if  the  gate  voltage  is  such 
that  there  is  a  charge  equivalent  to  one  half  electron  on  the  gate  capac¬ 
itance,  Cg  Vg  =  ( q0  +  })  e  mod  e,  tunnel  events  are  energetically  allowed 
(Fig.  4b).  A  cascade  of  transitions  between  charge  states  occurs,  and 
current  flows  between  source  and  drain. 

Note  that  in  practice  the  current  cannot  be  turned  off  completely 
in  the  Coulomb-blockaded  state.  Higher-order  processes  in  which 
several  tunnel  events  occur  simultaneously  (co-tunnelling  processes), 
and  whose  relative  importance  scale  as  powers  of  RK/RT,  will  induce  a 
small  leakage  current  not  described  by  the  orthodox  theory.  However, 
as  long  as  RK  RT,  the  SET  behaves  as  a  charge  amplifier  as  the  pres¬ 
ence  or  absence  of  a  fraction  of  an  electron  on  the  gate  capacitance  can 
control  an  easily  measurable  current  (the  order  of  magnitude  is 
around  109  e/s  in  practical  cases)  (Fig.  5).  This  gain  is  not  the  only 
factor  determining  the  sensitivity  of  the  device.  One  must  take  into 
account  the  shot  noise  in  the  source-drain  current,  which  is  attribut¬ 
able  to  the  quantum  randomness  of  the  time  intervals  between  tunnel 
events  (we  work  in  a  regime  where  thermal  fluctuations  can  be 
neglected).  This  quantum  randomness,  corresponding  to  electrons 
having  to  ‘choose’  between  the  two  sides  of  the  barrier,  is  the  process 
that  ultimately  limits  the  sensitivity  of  the  device.  As  we  will  show 
below,  the  noise  characteristics  of  the  SET  can  be  calculated  exactly  in 
a  simple  regime,  which  is  rich  enough  to  yield  semi-quantitative 
understanding  of  the  quantum  limit  of  sensitivity. 


Returning  to  the  offset  gate  charge,  it  is  believed  that  its  value  is 
determined  in  part  by  differences  between  the  work  functions  of  the 
metal  of  the  island  and  that  of  the  other  electrodes.  Even  minute 
variations  in  the  work  functions  are  sufficient  to  cause  q0  to  fluctuate  by 
numbers  much  greater  than  unity,  which  is  the  case  observed  in  prac¬ 
tice.  Another  potential  source  of  fluctuation  in  offset  charge  is  charge 
motion  in  the  substrate  or  even  in  the  tunnel  barrier  oxide.  The  absence 
of  control  over  offset  charges  severely  hinders  any  application  of  SETs 
to  digital  electronics,  where  the  gate  thresholds  must  be  rigorously 
fixed  and  uniform.  However,  SETs  can  still  be  used  as  sensitive  ampli¬ 
fiers  in  the  audio-frequency  (a.f.)  or  radio-frequency  (r.f.)  domains,  as 
a  simple  additional  feedback  circuit  can  compensate  for  the  fluctua¬ 
tions  in  offset  charge  that  occur  mainly  at  lower  frequencies. 

Noise  characteristics  of  an  amplifier 

Before  we  examine  how  quantum  shot  noise  affects  the  ultimate 
performances  of  SETs,  it  is  useful  to  discuss  the  noise  properties  of  a 
general  amplifier. 

A  linear  amplifier  can  be  described  phenomenologically  as  a 
‘black  box’  with  two  input  leads  and  two  output  leads.  We  can  repre¬ 
sent  the  inside  of  the  black  box  by  effective  elementary  components 
which  accurately  describe  how  it  appears  to  the  outside  circuitry.  If 
we  limit  ourselves  to  the  simpler  case  where  the  amplifier  has  no 
feedback  (that  is,  input  current  is  independent  of  output  voltage), 
we  arrive  at  the  schematics  of  Fig.  6  (feedback  introduces  complica¬ 
tions  that  are  not  crucial  in  our  discussion) .  Three  elements  describe 
the  transformation  of  the  signal  by  the  amplifier:  an  ideal  voltage 
amplifier  with  infinite  input  impedance,  zero  output  impedance 
and  a  voltage  gain  G( to),  and  the  input  and  output  impedances, 
2in(  a>)  and  Zout(  a>) .  In  these  parameters,  the  argument  co  denotes  the 
signal  angular  frequency. 

In  addition,  the  random  fluctuations  due  to  the  amplifier  are 
described  by  two  noise  sources  with  very  different  roles.  The  voltage 
noise  source  VN  describes  the  output  noise,  that  is,  the  noise  added  by 
the  amplifier  to  the  output  signal,  referred  to  the  input.  The  current 
noise  IN,  on  the  other  hand,  describes  the  back-action  of  the  amplifier 
on  the  circuit  at  the  input.  The  voltage  and  current  noise  are  assumed 
to  be  gaussian,  and  are  characterized  by  the  spectral  densities  Sy(co) 
and  SI(tt)),  which  are  the  Fourier  transforms  of  the  autocorrelation 
functions  of  VN  and  IN,  respectively.  We  neglect  here  the  correlations 
Sy^o)),  which  are  not  essential  for  our  discussion.  These  spectral 
densities,  together  with  the  input  impedance,  determine  the  ultimate 
resolving  power  of  the  amplifier  for  small  signals. 

It  is  useful,  in  the  discussions  of  this  resolving  power,  to  introduce 
several  combinations  of  these  quantities.  The  first  one  is  the  charge 
sensitivity  SQ(co)  =  V Sv(&/)/|ft>.Zto(&))|,  which  we  have  already  dis¬ 
cussed.  The  second  quantity  is  the  energy  sensitivity  e(co)  = 
SQ2(co)/2Cj]1  where  =  Re{  1  /[icoZ^co)]}.  Assuming  that  the  ampli¬ 
fier  is  driven  by  a  voltage  source  with  strictly  zero  internal  impedance 
(and  hence  always  in  the  classical  regime,  as  the  zero-point  voltage 
fluctuations  are  proportional  to  the  real  part  of  the  impedance),  this 
quantity  tells  us  the  smallest  amount  of  energy  SE  =  e/  r  fed  by  the 
source  into  the  input  circuit  of  the  amplifier  which  will  give  an  output 
signal  that  is  distinguishable  from  zero  when  we  accumulate  it  during 
a  time  r.  Even  though  e  has  the  dimension  of  an  action  and  is  conve¬ 
niently  measured  in  units  of  h,  quantum  mechanics  does  not  impose 
any  restriction  on  the  ratio  elfl,  which  can  in  principle  tend  to  zero23. 
Nevertheless,  in  systems  studied  so  far,  a  value  of  elh  of  order  unity 
usually  means  that  the  contribution  of  thermal  fluctuations  to  the 
noise  of  the  amplifier  have  been  suppressed  down  below  those  of 
quantum  fluctuations. 

Although  the  energy  sensitivity  seems  initially  to  be  a  very 
powerful  concept,  it  neglects  the  back-action  of  the  amplifier  and  is 
insufficient  to  determine  how  well  the  amplifier  performs  when  we 
use  it  to  measure  a  quantum  signal  coming  from  a  system  with  a  finite 
source  impedance  Zs.  The  current  noise,  which  induces  back-action 
voltage  fluctuations  across  the  source  impedance,  and  hence  on  the 
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Figure  5  Variation  at  T=  0  of  the  source-drain  current  in  a  SET  as  a  function  of  the 
voltage  between  drain  and  source  and  the  voltage  between  gate  and  source.  We  have 
assumed  that  the  two  junctions  have  identical  parameters.  The  sum  of  the  junction 
tunnel  resistances  is  twice  the  resistance  quantum  /7/e2  and  a  small  amount  of  co¬ 
tunnelling  rounds  the  peaks  associated  with  the  lifting  of  the  Coulomb  blockade  (gate 
charge  corresponding  to  half-integers  multiplied  by  the  charge  quantum).  The 
sensitivity  of  the  SET  is  based  on  the  rapid  variation  of  the  source-drain  current  as  the 
gate  charge  varies  by  a  fraction  of  one  electron. 


signal  going  into  the  amplifier,  must  be  considered  together  with  the 
voltage  noise. 

We  therefore  introduce  the  noise  impedance  Z N  =  V  Sy/S,  and  the 
noise  energy  £N  =  VS^Sj.  They  have  the  following  meaning:  suppos¬ 
ing  that  1/Zin(<u)  =  0,  the  amplifier  will  add  a  minimal  amount  of 
noise  to  the  signal  coming  from  the  source  when  Zs  =  Z^  (noise 
matching) .  This  minimal  noise  has  a  power  per  unit  bandwidth  given 
by  £n  (this  quantity  has  the  dimension  of  an  energy,  hence  the  name 
noise  energy).  Quantum  mechanics  places  a  strict  restriction  on  the 
noise  energy.  No  amplifier  can  have  an  smaller  than  hco/2,  half  a 
photon  energy  at  the  signal  frequency  co  (refs  23-25).  This  funda¬ 
mental  limitation  is  a  form  of  Heisenberg’s  uncertainty  principle  (see 
Table  1 ) .  It  is  worth  mentioning  that  in  the  classical  regime,  that  is,  the 
regime  where  the  contribution  of  quantum  fluctuations  to  the  noise 
of  the  amplifier  is  negligible,  the  noise  temperature  TN  =  E^l  kB  is  often 
used  in  place  of  the  noise  energy.  The  meaning  of  TN  corresponds  to  a 
well-defined  procedure:  imagine  one  connects  a  resistor  with  value 
Z^  at  the  input  of  the  amplifier.  The  temperature  to  which  one  must 
heat  the  resistor  to  double  the  noise  measured  at  the  output  of  the 
amplifier  is  precisely  TN. 

If  we  now  apply  these  concepts  to  FET  amplifiers,  we  find  that  the 
best  performance  in  the  r.f.  domain  0.1-10  GHz  is  obtained  with 
heterostructure  FETs  cooled  to  4  K  (ref.  26).  At  500  MHz,  their  noise 
temperature  is  around  2  K  (the  corresponding  thermal  energy  is 
equivalent  to  the  energy  of  40  photons  at  1  GHz),  their  noise  imped¬ 
ance  is  about  50  Cl  and  their  energy  sensitivity  is  of  the  order  102  h. 
They  are  thus  far  from  the  quantum  limit  (at  higher  frequencies  the 
minimum  noise  energy  improves  and  can  be  equivalent  to  only  ten 
photons).  As  a  charge-sensing  device,  their  best  performance  is 
around  10-2  e/VHz  (ref.  8). 

Although  theoretically  the  heterostructure  FET  could  reach 
the  quantum  limit  of  sensitivity  in  the  form  of  a  ballistic,  two- 
dimensional  electron  gas  quantum-point  contact27,  it  is  difficult  in 


Figure  6  Effective  circuit  elements  describing  the  properties  of  a  linear  voltage 
amplifier  with  no  feedback.  The  triangle  represents  an  ideal  noiseless  voltage  amplifier 
with  infinite  input  impedance  and  zero  output  impedance.  Its  gain  is  denoted  by  G.  The 
boxes  marked  2j„  and  ^correspond  to  the  input  and  output  impedances,  respectively. 
The  noise  sources  l/N  and  /N  represent  the  voltage  noise  (noise  added  by  the  amplifier 
at  the  output,  but  referred  to  the  input)  and  the  current  noise  (back-action  noise  of  the 
amplifier  sent  to  the  circuit  feeding  the  input).  These  elements  are  eventually 
dependent  on  the  signal  frequency  o> 


this  device  to  achieve  effective  electrostatic  coupling  of  the  gate 
electrode  to  the  channel  without  affecting  its  mobility  and  inducing 
parasitic  input  capacitance.  In  the  SET,  by  contrast,  there  is  more  flex¬ 
ibility  in  the  design  of  the  gate  capacitance.  The  fabrication  of  the 
latter  is  to  a  large  extent  disconnected  from  that  of  the  tunnel  barriers, 
whose  quality  is  analogous  to  the  mobility  of  the  channel  in  the  FET. 
We  will  now  examine  how  the  operating  principle  of  the  SET  brings 
us  close  to  the  quantum  limit. 

Ultimate  sensitivity  of  SETs 

In  the  framework  of  the  orthodox  theory,  the  tunnel  events  constitute 
a  generalized  Poisson  process  with  correlations.  It  is  possible  to  calcu¬ 
late  analytically  the  noise  characteristics  of  the  SET  at  temperatures 
such  that  kBT  <§  eVds  (refs  7,  28).  The  following  expressions  are 
obtained  for  the  voltage  and  current  noise,  in  the  simple  case  where 
the  two  junctions  have  identical  parameters: 


Sv(cu)  = 


(1  —  a2)(l  +  a2) 
8a2 


 ( 1  —  a2)  eyHs 

4 


In  what  follows,  we  neglect  the  influence  of  the  correlation 


Svi(w)  = 


(1  -a2) 


eVA. 


between  voltage  and  current  noise,  whose  effect  would  be  to  reduce 
the  noise  energy  if  an  appropriate  complex  source  impedance  could 
be  presented  to  the  SET. 

In  these  noise  expressions  the  Coulomb  blockade  parameter  a  = 
(2CgVg—  e)/C2Vds,  which  is  limited  here  to  the  range  0  <  a  < 
1  —  Max(_RK/77Gs,  eVdJkBT),  fixes  the  relative  values  of  the  gate  and 


Table  1  Constraints  on  sensitivity  and  back-action  imposed  by  quantum  mechanics 

System 

Sensitivity 

Back-action 

Limiting  relation 

Limited  quantity 

Heisenberg  microscope 

Ax 

A P 

AxAp  5=  h/2 

Action 

Electronic  amplifier 

Sv 

S, 

(SvS,)1'2  3*  fao/2 

Noise  energy 

Qubit  read-out 

7m 

G 

GG  9  1/2 

Information 

Table  1  shows  different  forms  of  Heisenberg’s  uncertainty  relation  linking  the  sensitivity  of  a  measurement  of  a  given  physical  quantity  and  the  simultaneous  back-action  on  the  conjugate  quantity.  The 
precision  with  which  the  position  of  an  object  would  be  measured  with  a  photon  is  related  to  the  momentum  transferred  to  this  object  by  the  radiation  pressure  of  the  photon  (row  1).  Similarly,  the 
voltage  sensitivity  with  which  an  amplifier  measures  the  circuit  at  its  input  is  related  to  the  back-action  current  noise  that  this  amplifier  produces  in  the  circuit  (row  2).  Finally,  the  time  needed  to  acquire 
the  value  of  a  qubit  is  related  to  the  rate  of  dephasing  of  the  qubit  induced  by  the  back-action  of  the  read-out  (row  3).  Quantities  characterizing  sensitivity  and  back-action,  respectively,  are  given  in 
columns  2  and  4. 
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drain-source  voltage.  For  1  —  RK/(^ rRf)  <  a,  co-tunnelling  processes 
dominate  over  the  single  tunnelling  processes  we  consider  here,  and 
for  1  —  eVislkBT  <  a  the  influence  of  thermal  fluctuations  would  to 
have  be  taken  into  account  in  the  above  expressions.  The  resistance 

=  RT1  +  RT1  is  the  sum  of  the  two  junction  resistances.  The  above 
expressions  are  valid  for  source-drain  voltages  below  the  Coulomb 
gap  VjsCj/e  <  1  and  at  frequencies  co  that  are  low  on  the  tunnel  rate 
scale  (eRj/Vj,)-1.  We  have  taken  the  mean  offset  charge  to  be  zero 
because  it  appears  only  as  a  shift  in  Vg,  but  the  typical  value  of  the 
expected  Ilf  offset  charge  fluctuations  make  our  expressions  valid 
only  above  the  crossover  to  the  intrinsic  device  shot  noise  at  a  few  105 
Hz.  Note  that  at  this  level  of  approximation  the  input  impedance  is 
simply  a  capacitance  resulting  from  the  series  combination  of  the  gate 
and  the  sum  Ct  =  Cn  +  C,2  of  the  two  junction  capacitances:  = 

Czl[io)Cg(Cn  +  CJ2)].  A  dissipative  part  in  the  input  impedance 
appears  only  at  higher  orders  in  the  dimensionless  signal  frequency 
(ewiyvj-1. 

It  is  straightforward  to  go  from  these  expressions  to  the  charge 
sensitivity,  energy  sensitivity,  noise  energy  and  noise  impedance: 


«CX®)  = 


e(w)  = 


Exit o)  = 


C,  / 
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Rk  \elCs  \C, 
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2  a 


2  Rr 


ha> 


Zs(<0)  =  . 


1  +  of  /  Q,V 


8a"  l  CJ  eu> 
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C„E, 
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We  can  find  a  conservative  upper  bound  for  the  optimal  value  of 
the  quantities  that  characterize  the  sensitivity  of  the  SET  by  minimiz¬ 
ing  the  above  expressions  over  the  range  of  parameter  values  that 
correspond  to  our  hypotheses.  We  therefore  take  a  =  1  —  2RK/3R-Zand 
Rt  =  2 Rk,  values  at  which  co-tunnelling  remains  marginal.  We  also 
take  the  source-drain  voltage  Vds  =  el(2Cf)  to  be  able  to  neglect  ther¬ 
mal  fluctuations  in  practical  situations.  In  addition,  we  take  Cj  =  1  fF 
as  a  compromise  between  reaching  attainable  temperatures,  keeping 
the  heating  of  the  island  by  the  drain-source  current  at  a  reasonable 
level  and  achieving  an  acceptable  output  signal  (in  contrast  with  the 
previous  parameter  choices,  this  last  value  is  dictated  essentially  by 
the  current  technology  issues  and  not  by  the  validity  of  our  expres¬ 
sions)  .  We  arrive  at  the  optimal  values 

SQopt—  1-7  X  10-6e/VHz 
eopt£0.7  h(CtICs) 

Tjsjopt  ~  2.2ho) 

[Gx  EN]opt  =  0.14eVdsS  75 mKXfcB 

We  have  left  the  internal  coupling  ratio  (Cj/Cg)  of  the  SET  in  the 
right-hand  side  of  the  expression  for  eopt  as  this  parameter  can  easily 
be  modified  by  changing  the  lithography  of  the  device.  This  ratio 
determines  how  much  of  the  energy  fed  by  the  gate  voltage  into  the 
SET  is  used  to  charge  the  island,  which  is  the  active  element,  rather 
than  the  gate  capacitance. 

The  value  for  EM,  shows  that  the  SET  operates  in  the  vicinity  of 
the  quantum  limit.  The  product  [GxfN]op,  gives  the  SET  output 
noise  which  would  ideally  be  higher  than  the  added  noise  of  the 
following  amplifier,  that  is,  a  cryogenic  heterostructure  FET.  Today, 
this  FET  is  the  limiting  factor  in  the  system  noise,  and  the  best  experi¬ 
mental  upper  bound  so  far  for  the  energy  sensitivity  is  40fz(Cj/Cg) 


Figure  7  The  single-Cooper-pair  box.  a,  Circuit  diagram  of  the  single-Cooper-pair  box. 
This  simple  system  consists  of  one  superconducting  tunnel  junction  (box  with  cross)  in 
series  with  a  voltage  source  L/and  a  gate  capacitance  Cg.  When  the  charging  energy  of 
the  island  is  less  than  the  superconducting  gap,  all  electrons  are  paired  and  the  number 
n  of  excess  electrons  is  even.  Cooper  pairs  can  tunnel  reversibly  through  the  tunnel 
junction.  Charge  states  differing  by  one  Cooper  pair  are  coupled  by  the  Josephson 
energy  £,.  b,  Energy  of  the  two  lower  quantum  states  of  the  Cooper-pair  box  (full  line)  as 
a  function  of  U.  The  dashed  lines  represent  the  electrostatic  energy  of  the  box  including 
the  work  done  by  the  voltage  source.  At  the  avoided  crossings,  the  two  charge  states 
are  mixed  and  constitute  a  solid-state  implementation  of  a  qubit, 
c,  Quantum  mechanical  average  of  the  number  of  electrons  in  the  ground  state  of  the 
box,  as  a  function  of  U. 


in  the  r.f.  domain4.  The  improvement  of  the  combination  of  the  SET 
with  a  FET,  or  with  a  superconducting  quantum  interference  device 
(SQUID)  amplifier  (see  below) ,  as  well  as  the  verification  of  the  above 
predictions  for  the  noise,  is  a  topic  for  future  research. 

Other  quantum-limited,  solid-state  r.f.  amplifiers 

According  to  the  above  analysis,  the  SET  approaches  but  does  not 
quite  reach  the  quantum  limit.  The  back-action  noise  due  to  the 
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Preparation  Free  evolution  Measurement 


fW 


P2=(1-exp(fw/f0)) 


Figure  8  Measuring  the  state  of  a  quantum  system,  a,  Schematic  circuit  showing  how 
a  qubit  implemented  by  a  Cooper-pair  box  could  be  read-out  by  an  r.f.-SET.  b,  Time 
evolution  of  signals  in  the  circuit  of  a.  While  the  read-out  r.f.-SET  is  turned  off,  a  first 
pulse  prepares  the  qubit  in  a  quantum  superposition  of  ground  and  excited  state  at  the 
crossing  point  CJU=e.  After  a  waiting  period,  a  second  pulse  is  applied  which,  in 
absence  of  decoherence,  would  make  the  qubit  return  to  the  ground  state.  Finally,  after 
these  manipulations,  the  read-out  transistor  is  switched  on  while  the  box  gate  charge  is 
moved  away  from  Csil=e.  If  decoherence  occurs,  the  qubit  will  be  left  in  the  n = 2 
charge  state  with  a  non-zero  probability  p2.  Given  its  sensitivity  and  low  back-action 
noise,  the  r.f.-SET  should  be  able  to  detect  this  event  with  a  signal-to-noise  ratio  better 
than  unity.  Ct,  and  L„  refer  respectively  to  the  capacitance  and  inductance  of  the 
impedance  transformer.  The  times  fw  and  are  the  waiting  time  between  the  two 
pulses  and  the  intrinsic  decoherence  time  of  the  box,  respectively. 

randomness  of  tunnel  events  is  dominated  by  phase-scrambling 
processes  inside  the  island.  But  in  the  co-tunnelling  regime, 
back-action  noise  is  dominated  by  fluctuations  of  the  island  voltage 
associated  with  the  measurement  of  the  gate  voltage  itself  and  should 
therefore  be  more  efficient.  More  work,  experimental  and  theoreti¬ 
cal,  is  needed  to  quantify  this  conjecture. 

The  superconducting  SET  can  measure  the  charge  on  the  gate 
with  fully  coherent  carriers,  the  Cooper  pairs.  Zorin  has  analysed 
theoretically  a  version  of  this  device  shunted  by  a  resistor  and  found 
that  the  noise  energy  could  in  this  case  exactly  reach  the  quantum 
limit29.  However,  the  mode  of  operation  seems  to  involve  a  more 
precise  tuning  of  the  device  parameters  than  for  the  SET. 

So  far  we  have  considered  voltage  amplifiers,  that  is,  amplifiers 
with  an  input  impedance  that  is  large  compared  with  the  resistance 
quantum  RK.  For  several  years  already,  SQUID  devices  have  played 
the  role  of  the  SET  in  the  realm  of  low-impedance  amplifiers.  There 
is  in  fact  a  duality  relationship  between  the  d.c. -SQUID  and  the 
SET1.  Whereas  the  SET  is  based  on  charge  quantization  of  a  metallic 
island  surrounded  by  an  insulator,  the  d.c.-SQUID  is  based  on  the 
quantization  of  flux  in  a  superconducting  loop.  The  sub-electron 
sensitivity  of  the  SET  has  its  analogue  in  the  sub-flux  quantum  sensi¬ 
tivity  of  the  SQUID.  Although  it  has  been  known  for  many  years  that 
a  SQUID  could  operate  at  the  quantum  limit,  this  has  been  achieved 
only  recently  in  the  r.f.  domain.  Andre  et  al.  have  shown  that  a 
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SQUID  with  a  microstrip  input  line  could  achieve  A:BTN  —  5hco  at 
438  MHz  (ref.  30). 

Finally,  we  should  mention  a  special  class  of  linear  amplifiers: 
the  mixers  based  on  photon-assisted  tunnelling  in 
superconductor-insulator-superconductor  (SIS)  junctions. 
They  convert  a  signal  at  several  hundred  GHz  into  a  signal  at  several 
GHz.  Even  though  the  absolute  power  of  the  output  signal  is 
weaker  than  the  power  of  the  input  signal,  these  devices  have  a 
large  ‘photon  number’  gain.  Futhermore,  their  noise  closely 
approaches  the  quantum  limit:  EN  =  0.6f7<w(ref.  31).  Closely  related  to 
SIS  mixers  are  Josephson  parametric  amplifiers,  which  have  been 
operated  at  the  quantum  limit32.  However,  these  devices  do  not 
have  the  advantage  of  the  SET  and  the  SQUID  to  amplify  at  the 
same  time  both  a.f.  and  r.f.  signals,  a  useful  feature  for  tuning  and 
trouble-shooting. 

Measuring  the  state  of  a  two-level  system 

In  the  past  few  years  there  has  been  much  interest  in  the  possibility  of 
realizing  a  quantum  computer6.  Although  the  more  advanced 
experiments  in  this  field  are  taking  place  in  quantum  optics 
systems33’34,  several  implementations  of  quantum  bits  and  quantum 
gates  in  solid-state  systems  have  been  proposed35-38.  We  focus  here  on 
a  charge  qubit  which  is  now  well  understood  experimentally12,13  and 
theoretically39,40:  the  Cooper-pair  box  (Fig.  7).  It  consists  of  a 
superconducting  island  with  Coulomb  energy  Ec  connected  to  a 
superconducting  charge  reservoir  by  means  of  a  Josephson  tunnel 
junction  with  Josephson  coupling  energy  £,.  The  island  is  influenced 
electrostatically  through  a  gate  capacitor  Cg  connected  to  voltage 
source  (Fig.  7a).  If  the  conditions  kBT  <4  Ej  Ec  <  A  are  realized, 
where  A  is  the  energy  gap  of  the  superconductor,  then  the  island 
will  have  only  an  even  number  n  of  excess  electrons.  This  only  degree 
of  freedom  will  be  a  good  quantum  number,  except  at  the  avoided 
crossings  of  the  charge  levels  (Fig.  7b, c).  A  recent  experiment  has 
shown  that  such  a  solid-state  qubit  was  able  to  display  several  Rabi 
oscillations  when  stimulated  by  a  microwave  pulse13.  However,  the 
measurement  scheme  in  this  experiment  involved  a  probe  tunnel 
junction  which  observed  the  qubit  continuously,  but  weakly.  Such  a 
continuous  observation  dephases  the  qubit  and  prevents  a 
measurement  of  its  intrinsic  decoherence  time.  It  is  therefore 
desirable  to  read-out  the  qubit  only  at  the  end  of  the  coherent 
manipulation  and  evolution  period  (Fig.  8).  An  obvious 
candidate  for  the  read-out  is  the  SET.  In  a  version  with  r.f.  bias4,  the 
transistor  can  be  turned  on  and  off  very  quickly,  the  turn-on 
time  being  set  by  the  damping  time  of  the  tank  circuit.  This  latter 
is  of  the  order  of  several  tens  of  nanoseconds,  a  time 
supposedly  shorter  than  the  intrinsic  decoherence  time  of  the  box, 
estimated  to  be  longer  than  1  ps.  The  question  therefore  arises  as  to 
whether  the  SET  is  sufficiently  sensitive  to  detect  the  state  of  the 
box  with  a  signal-to-noise  ratio  allowing  the  study  of  decoherence 
mechanisms. 

Measurement  time  and  dephasing  rate 

An  important  concept  is  the  measurement  time  Tm  of  the  read-out 
SET14,  which  is  defined  as  the  time  needed  for  discriminating 
between  the  two  charge  states  of  the  box  differing  by  a  Cooper  pair, 
assumed  to  be  good  eigenstates,  with  a  signal-to-noise  ratio  of  1. 
Island  box  charge  is  a  good  quantum  number  when  EiIEA  <§  1,  where 
the  electrostatic  energy  Ed  =  2e2(l  —  CgU/e)/Ctot  of  the  box  is 
controlled  by  the  voltage  U  (here  Cg  and  Ctot  refer  to  the  box  gate 
capacitance  and  total  island  capacitance).  In  practice,  the  charge 
measurement  is  made  just  after  U  is  varied  away  from  the  crossing 
point  (Fig.  8b).  The  circuit  of  Fig.  8a  shows  that  Tm  is  closely  related  to 
the  voltage  noise  of  the  SET,  with  Tm  =  4Sv/(2e/Ctot)2.  We  thus  find 
that  Tm  =  (5Q/e)2(Ctot/Cta)2,  where  CJCtot  is  the  charge-coupling 
coefficient. 

It  is  interesting  to  note  that  the  back-action  current  noise 
dephases  the  charge  states  relative  to  one  another  only  if  charge  is  a 
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good  quantum  number.  The  general  expression  for  the  dephasing 
rate  is 


hi  £j  +£el 


Sv°x(&>=0) 


where  Sy™  (<w  )  =  S^co)/ Ctotco)1 2 3 4  is  the  fluctuation  spectral  density  of 
the  box  island  voltage  induced  by  the  SET  current  noise.  In  the  regime 
£,/-Eel  1,  we  find  that  the  product  TmQ,  =  2 [E^lcol/hco]^  involves 

only  the  noise  energy  and  is  of  order  unity.  Once  again,  this  is  another 
close  approach  to  the  quantum  limit  as  the  minimum  value  for  Tm  Fr() 

is|  (Table  !)• 

So  in  principle  it  seems  that  the  SET  could  acquire  charge  with  any 
given  large  signal-to-noise  ratio  by  measuring  the  box  for  a  long 
enough  time.  However,  any  coupling  between  charge  states  will 
induce  transitions  and  will  corrupt  the  measurement.  The  transition 
rate  between  charge  states  due  to  the  current  noise  is  given  by 


r,= 


hi  E2,+Ee i 


Sv°x(ft)  =  /2) 


where  Cl  =  V E2j  +  E2Jh  is  the  transition  frequency  between  charge 
states. 

We  thus  find  that  the  signal-to-noise  ratio  including  the  effect  of 
transitions  between  charge  states  is 


s/n=  (r1rj_1/2 


1  /  2E;  \~m 

( ENlhco)m=0\E)+E2el ) 


In  practice,  a  signal-to-noise  ratio  significantly  greater  than  1  seems 
possible  with  an  optimized  r.f.-SET  in  a  single  charge  measurement, 
even  when  taking  into  account  the  50%  reduction  in  sensitivity  of  the 
SET  in  the  r.f.-bias  mode41. 

Towards  single-photon  sensitivity  in  the  r.f.  domain 

To  summarize,  the  SET  transistor  is  a  charge  amplifying  device  whose 
sensitivity  in  the  r.f.  domain  is  limited  in  principle  only  by  a  well 
understood  process,  quantum  shot  noise.  This  property  displays  a 
marked  contrast  with  a  FET.  It  should  be  possible  in  the  near  future  to 
show  experimentally  that  the  noise  energy  of  the  SET  approaches  the 
quantum  limit  within  a  factor  of  order  unity.  Although  the  SET  will 
not  be  100%  efficient  in  acquiring  charge  information  with  only  the 
minimal  back-action  noise  required  by  quantum  mechanics,  it 
should  still  be  able  to  read-out  a  charge  qubit  in  a  single-shot 
measurement.  In  more  sophisticated  devices  using  Cooper-pair 
tunnelling  or  co-tunnelling  processes,  the  quantum  limit  could  be 
approached  even  more  closely.  Furthermore,  by  transposing  in  the 
r.f.  domain  the  manipulations  of  the  quantum  signal  that  are  now 
performed  routinely  in  experiments  in  cavity  quantum  electrody¬ 
namics42,  one  could  use  these  ultimate  amplifiers  for  measuring 
signals  consisting  of  a  single  photon,  without  even  destroying 
the  photon.  □ 
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Abstract 

We  demonstrate  that  a  continuously  measured  microelectronic  circuit,  the  Cooper-pair  box  mea¬ 
sured  by  a  radio-frequency  single-electron  transistor  (RF-SET),  approximates  a  quantum  two-level 
system.  We  extract  the  Hamiltonian  of  the  circuit  through  resonant  spectroscopy  and  measure 
the  excited-state  lifetime.  The  lifetime  is  more  than  105  times  longer  than  the  inverse  transition 
frequency  of  the  two-level  system,  even  though  the  measurement  is  active.  This  lifetime  is  also 
comparable  to  an  estimate  of  the  known  upper  limit,  set  by  spontaneous  emission,  for  this  circuit. 

PACS  numbers:  74.40+k,  85.25Na,  85.35Gv 


1 


Recently,  microelectronic  circuits  have  been  coaxed  into  behaving  as  quantum  two-level 
systems  (TLS)  [1-5].  The  TLS  behavior  of  circuits  is  revolutionary  because  it  demonstrates 
the  quantum  behavior  of  a  macroscopic  degree  of  freedom  composed  of  many  microscopic 
degrees  of  freedom.  Quantum  coherence  was  believed  to  be  fragile  in  electrical  circuits 
both  because  it  required  the  suppression  of  the  dynamics  of  the  microscopic  elements  in  a 
condensed  matter  system,  and  because  the  quantum  oscillations  of  an  electric  or  magnetic 
degree  of  freedom  would  efficiently  radiate  energy  into  the  electromagnetic  environment. 
Discussed  in  terms  of  the  Bloch  equations  [6],  familiar  from  nuclear  magnetic  resonance,  a 
TLS  in  a  coherent  superposition  of  states  has  characteristic  times  T2  to  become  an  incoherent 
mixture  and  Ti  to  relax  back  to  its  ground  state. 

In  this  paper,  we  observe  that  a  microelectronic  circuit,  the  Cooper-pair  box,  may  be 
measured  continuously  while  still  behaving  approximately  as  a  two-level  system.  The  box 
is  integrated  with  a  measurement  apparatus,  the  RF-SET,  which  we  operate  as  weak,  con¬ 
tinuous  measurement  of  the  box’s  state.  Under  these  conditions  we  are  able  to  determine 
the  parameters  that  appear  in  the  box’s  Hamiltonian,  make  a  worst  case  estimate  T2*  of  the 
decoherence  time  T2,  and  measure  the  excited-state  lifetime  Tf  of  the  two-level  system.  We 
determine  the  parameters  in  the  Hamiltonian  through  a  kind  of  spectroscopy  where  we  ob¬ 
serve  a  resonant  change  in  the  box’s  state  when  its  transition  frequency  matches  a  multiple 
of  the  frequency  of  an  applied  oscillatory  excitation.  From  the  width  in  frequency  of  these 
resonances  we  can  find  T2  [7].  We  stimulate  the  box  into  its  excited  state  and  measure 
T\  directly  by  exploiting  the  large  measurement  bandwidth  of  the  RF-SET  to  resolve  in 
time  the  circuit’s  decay  to  its  ground  state.  Most  remarkably,  the  value  of  T\  that  we  find 
while  continuously  measuring  the  state  of  the  box  is  comparable  to  estimates  of  the  excited- 
state  lifetime  limited  by  the  quantum  fluctuations  of  the  electromagnetic  environment.  This 
demonstrates  that  the  Cooper-pair  box,  when  embedded  in  a  circuit  for  control  and  mea¬ 
surement,  remains  well  decoupled  from  other  sources  of  dissipation.  Based  on  the  observed 
noise  in  the  readout  and  the  lifetime  we  conclude  that  RF-SET  is  a  promising  qubit  readout 
because  a  ’single-shot’  measurement,  where  the  box  is  observed  in  its  excited  state  before  it 
has  relaxed  into  its  ground  state,  is  possible. 

The  Cooper-pair  box  is  a  microelectronic  circuit  composed  of  an  isolated  superconducting 
island,  attached  to  a  superconducting  lead  through  a  thin  insulating  layer  across  which 
Cooper-pairs  can  tunnel.  An  additional  lead,  called  the  gate  lead,  lies  near  the  island  and 
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FIG.  1:  (a)  An  SEM  micrograph  of  the  Cooper-pair  box  and  SET  electrometer.  The  device  is 
made  from  an  evaporated  aluminum  film  (light  gray  regions)  on  an  insulating  Si02  substrate  (dark 
gray  regions)  by  the  technique  of  double  angle  evaporation  [8],  which  gives  the  double  image.  The 
aluminum  has  BCS  gap  A /ks  =  2.4  K.  (b)  A  circuit  diagram  of  the  box  and  RF-SET  electrometer 
showing:  the  voltage  Vg  and  magnetic  flux  $  which  control  the  box’s  Hamiltonian,  the  quantities 
VgC  and  u>  which  set  the  amplitude  and  frequency  of  the  microwave  excitation,  the  voltage  Vge  and 
the  dc  and  500  MHz  oscillatory  (RFin)  bias  which  determine  the  electrometer’s  operating  point, 
and  the  capacitance  Cc  that  couples  charge  between  box  and  electrometer.  The  charge  on  the  box 
is  inferred  from  variation  in  the  amount  of  applied  RF  power  that  is  reflected  (RFout)  from  the 
SET  electrometer,  which  is  a  sensitive  function  of  SET’s  conductance  [9].  The  tunnel  junctions 
(crosses  in  boxes)  are  characterized  by  a  junction  resistance  Rj  and  capacitance  Cj,  which  control 
the  box’s  Hamiltonian  through  CV  =  Cc  +  2 Cj  +  Cg  and  Rc,  =  Rj / 2  (see  text). 

changes  the  electrostatic  potential  of  the  island  with  the  application  of  a  voltage  Vg  to  the 
gate  lead  through  the  gate  capacitance  Cg  [Fig.  1(a)].  The  island’s  total  capacitance  CV; 
is  small  enough  that  the  addition  of  a  single  Cooper-pair  to  the  island  requires  a  large 
electrostatic  energy,  leading  to  suppressed  fluctuations  of  charge  on  the  island.  Because 
the  island  and  the  lead  are  superconducting,  all  of  the  electrons  form  Cooper-pairs  and 
participate  in  the  macroscopic  quantum  ground  state  of  the  island.  The  only  degree  of 
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freedom  is  the  number  of  pairs  n  on  the  island.  Because  of  the  large  charging  energy, 
we  need  only  consider  two  states,  a  state  |0)  with  no  excess  Cooper-pairs  (n  =  0),  and  a 
state  |1)  with  one  excess  Cooper-pair  {n  =  1),  as  reckoned  from  electrical  neutrality.  The 
Hamiltonian  of  the  Cooper-pair  box  circuit  is 

H  =  -2EC(1  -  2 ng)crz  -  (1) 

where  crz  and  crx  are  the  Pauli  spin  matrices  and  ng  is  total  polarization  charge  applied 
to  the  gate  electrode,  ng  =  CgVg/2e  —  n0ff,  in  units  of  a  Cooper-pair’s  charge  [10,  11]. 
The  offset  charge  n0ff  accounts  for  the  uncontrolled  potential  arising  from  charges  nearby 
the  box  island.  The  charging  energy,  Eq  =  e2/2 Cs,  is  the  electrostatic  energy  required 
to  add  one  electron  to  the  island  and  the  Josephson  energy,  E™^  =  hA/8e2R s,  is  the 
effective  tunnelling  matrix  element  for  Cooper-pairs  across  a  junction  with  resistance  Ry, 
in  a  superconductor  with  BCS  gap  A.  The  junction  is  in  fact  a  composite  of  two  parallel 
junctions  connected  to  form  a  loop  with  1  (/mi)2  area  (Fig.  1).  The  effective  Josephson 
energy  Ej  of  the  pair  of  junctions  is  then  tuneable  with  magnetic  flux  <3>  through  this 
loop,  as  Ej  =  Eyax  cos(7r<h/<l)0),  where  <h0  is  the  quantum  of  flux  (h/2e).  Equation  1  is  the 
Hamiltonian  of  a  quasi-spin  1/2  particle  in  a  fictitious  magnetic  field  that  can  be  decomposed 
into  two  orthogonal  fields.  The  2  component  of  this  fictitious  field  which  accounts  for  the 
box’s  electrostatic  energy,  Eei(Vg)  =  2EC(1  —  2 ng)  ,  is  tuned  with  Vg  and  the  x  component, 
which  accounts  for  the  Josephson  energy  E/(<f>)  =  Eyax  cos(7r<h/<J>o),  is  tuned  with  <J>  [11], 
The  box  is  an  artificial  two-level  system  and  both  of  the  terms  in  its  Hamiltonian  are  tuneable 
in  situ. 

In  the  box,  states  of  definite  numbers  of  Cooper  pairs  on  the  island  are  states  of  definite 
charge.  In  order  to  measure  the  charge  of  the  Cooper-pair  box,  we  fabricate  the  box  next 
to  a  radio-frequency  single-electron  transistor  (RF-SET)[8,  9],  an  exquisitely  sensitive  elec¬ 
trometer,  so  that  the  addition  of  a  Cooper-pair  to  the  box’s  island  causes  a  small  fraction 
(Cc/CC= 3.7%)  of  the  Cooper-pair’s  charge  to  appear  as  polarization  charge  on  the  capac¬ 
itor  Cc  that  couples  the  box  and  the  RF-SET  (Fig.  1).  The  electrometer  used  here  had  a 
sensitivity  of  4  x  1CT5  e/\/Hz  and  10  MHz  of  measurement  bandwidth.  Because  the  RF-SET 
measures  charge,  its  action  can  be  described  as  projecting  the  state  of  the  box  into  a  state  of 
definite  Cooper-pair  number.  In  the  formal  terms  of  Eq.  1,  it  measures  Qbox  =  (1  +  (erz))e 
where  Qbox  is  further  averaged  over  the  measurement  time. 
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FIG.  2:  (a)  The  ground  and  excited  state  energies  versus  ng  for  Eq.  1,  with  4 Ec  =  12 Ej  (solid 
line)  and  Ej  =  0  (dashed  lines).  Energy  eigenstates  asymptotically  approach  charge  states  (|1)  and 
|0»  far  from  ng  =  0.5.  (b)  Qbox  vs.  ng,  calculated  for  the  ground  state  (dotted  line),  excited  state 
(dashed  line),  and  measured  (solid  line)  with  35  GHz  microwaves  applied  to  the  box  gate.  The 
arrows  indicate  resonant  peaks.  Also  shown  is  Qbox  measured  with  no  microwaves  applied  (solid 
line),  with  the  y-axis  shifted  down  by  2.2  e.  (c)  Two  resonant  peaks  in  Qbox  vs.  ng  on  the  bottom 
axis  and  vs.  ujq\  on  the  top  axis,  with  u  =  38  GHz  and  where  the  larger  value  of  Vgc  (squares)  is 
twice  the  smaller  value  (triangles). 

We  perform  spectroscopy  by  applying  a  continuous  microwave  stimulus  to  the  gate  of  the 
Cooper-pair  box,  and  sweeping  ng  to  tune  the  parameters  of  the  TLS  and  find  the  resonance 
condition  (Fig.  2).  A  measurement  of  Qbox  vs.  ng  shows  that  the  box  does  not  remain  in 
its  ground  state  over  a  range  0.3  <  ng  <  0.7.  This  behavior  is  caused  by  backaction  [12,  13] 
generated  by  currents  flowing  through  RF-SET  [14].  We  proceed  by  studying  the  box  in  the 
range  of  ng  where  it  does  remain  in  its  ground  state. 

When  a  35  GHz  microwave  signal  is  applied  to  the  gate,  we  observe  clear  evidence  that 
the  box  is  a  coherent  two-level  system.  Resonant  peaks  appear  [Fig.  2(b)]  in  Qbox  that  are 
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sharp  and  symmetrically  spaced  about  ng  =  0.5.  The  two  features,  a  peak  for  ng  <  0.5  and 
a  dip  for  ng  >  0.5,  both  correspond  to  the  change  in  Qbox  when  the  box  spends  some  time  in 
the  excited  state.  Because  Qbox  is  an  average  of  thousands  of  repeated  measurements,  the 
peak  height  indicates  the  probability  of  finding  the  box  in  its  excited  state  [Fig  2(c)]. 

The  resonant  peaks  permit  a  spectroscopic  determination  of  Ec  and  E'} iax.  By  tuning 
rig  and  $  while  exciting  the  box  with  a  fixed  microwave  frequency,  we  find  good  agree¬ 
ment  between  the  locations  of  resonant  peaks  and  the  difference  E01(ng,  <f>)  =  hu0i  between 
ground-state  and  excited-state  energy  expected  from  Eq.  1.  An  independent  measurement 
of  Ec  [15]  demonstrates  that  these  peaks  occur  when  the  irradiating  frequency  to  is  half  ca0i, 
indicating  that  these  peaks  correspond  to  a  two-photon  transition  [16].  At  lower  frequencies 
and  for  single-photon  transitions,  the  peaks  would  appear  at  an  ng  for  which  the  box  does  not 
stay  in  the  ground  state  while  being  measured  and  are  therefore  not  visible.  We  End  a  single 
value  for  Ec  and  for  Ej mx  that  account  for  the  location  of  the  resonant  peaks  at  applied 
frequencies  between  32  and  38  GHz  giving  resonant  peaks  for  c^oi  between  64  and  76  GHz 
[Fig.  3(a)].  We  are  able  to  extract  the  parameters  of  the  Hamiltonian,  AEc/h  =  149.1  ±0.4 
GHz  and  Efax/ h  =  13.0  ±  0.2  GHz,  which  imply  CG  =  518  aF  and  =  12.1  kfl  Through 
spectroscopy  we  have  measured  the  parameters  of  an  electrical  circuit  that  could  not  have 
been  measured  with  transport  [Fig.  1(b)].  Because  these  measurements  were  made  at  a 
temperature  T  <  40  mK,  they  are  in  the  limit  ksT  <C  Ej  <  Ec- 

Consistent  with  the  behavior  of  a  TLS,  the  peaks  disappear  for  $  =  $0/2  when  Ej 
approaches  zero.  This  demonstrates  that  Ej  provides  the  coupling  between  the  charge  states 
[Fig.  3(b)].  An  oscillating  gate  voltage  with  amplitude  VgC  adds  a  term  to  the  Hamiltonian 
in  Eq.  1  which  is  (CgVgC /2e)  cos(ut)crZj  and  is  collincar  with  ground  state  of  the  quasi-spin 
described  by  Eq.  1  when  Ej  =  0.  The  microwave  excitation  therefore  applies  no  torque 
which  could  excite  the  quasi-spin  from  its  ground  state  [6]. 

The  width  of  the  resonant  peaks  we  observe  provides  a  worst-case  estimate  of  the  deco¬ 
herence  time  of  the  two-level  system.  We  express  the  width  of  a  resonance  5ng  as  a  width 
in  frequency  Anoi  =  ( l/h)(dE0i/dng)8ng .  In  the  absence  of  inhomogenous  broadening,  the 
half- width  at  half  maximum  inferred  for  zero  power  is  the  decoherence  rate  I/X2  of  a  TLS 
[6] .  From  5ca0i  measured  at  the  lowest  value  of  VgC  applied,  we  estimate  a  time  T%  of  about 
325  ps  [7].  The  resonant  peaks  have  a  Gaussian  shape,  and  nQff  drifts  an  amount  compara¬ 
ble  to  Sng  during  the  two  minutes  required  to  complete  a  measurement.  Both  observations 


6 


„  0.30 

-£  0.29 

|  0.28 
ns 

8  0.27 
^  0.26 
S.  0.25 


(<s>/o0) 


FIG.  3:  Resonant  spectroscopy  of  the  box  versus  the  two  control  parameters  of  the  Hamiltonian, 
Vg  and  $.  (a)  The  locations  of  resonant  peaks  (circles)  in  ng  and  4>,  for  to  =  32,  35,  and  38  GHz 
and  fits  (lines),  using  Eq.  1  for  =  2u ;  =  64,  70  and  76  GHz  to  find  single  value  of  Eq  and  of 
Ej-ax .  The  systematic  uncertainty  in  ng  is  represented  by  the  size  of  the  open  circle  symbols,  (b) 
The  height,  in  electrons,  of  a  76  GHz  resonant  peak  as  a  function  of  <h  (squares)  and  a  guide  to 
the  eye  (line). 

imply  that  the  width  of  the  peaks  expresses  not  the  intrinsic  loss  of  phase  coherence  dne 
to  coupling  the  TLS  to  the  environment,  but  rather  the  degree  to  which  an  ensemble  of 
measurements  are  not  identical,  due  to  the  well-known  1  //  noise  of  single-electron  devices 
[17]  .  This  T2*  is  a  worst-case  estimate  because  it  is  extracted  while  the  system  is  measured 
continuously  by  the  RF-SET  and  because  it  represents  an  ensemble  average  of  many  single 
measurements  that  require  about  two  minutes  to  complete.  Nevertheless,  T2  is  about  150 
times  longer  than  l/c^oi  [Fig.  2(c)]  and  is  similar  to  the  times  found  in  [18],  another  Cooper- 
pair  box  implementation,  as  well  as  [5]  a  SQUID  circuit.  Reference  [4]  demonstrates  that 
this  inhomogenous  broadening  may  be  overcome  by  operating  the  Cooper-box  at  ng  =  0.5 
where  Eq\  is  to  first-order  insensitive  to  fluctuations  in  naf / . 

In  order  to  measure  the  excited-state  lifetime  T±,  we  excite  the  box  and  then  measure  the 
time  required  to  relax  back  to  the  ground  state.  A  38  GHz  signal  is  continuously  applied  to 
the  gate  and  the  box  gate  is  tuned  to  ng  =  0.248  and  Ej  =  E™ax  so  that  the  microwaves 
resonantly  couple  the  ground  and  excited  state  through  a  two-photon  transition.  Abruptly, 


7 


FIG.  4:  A  determination  of  the  excited-state  lifetime  of  the  box.  Qbox  vs.  time  t,  (triangles), 
relative  to  t  =  0,  when  ng  is  shifted  from  =  0.248  to  0.171  in  30  ns,  with  38  GHz  microwaves 
applied.  The  shift  in  ng  brings  the  box  out  of  resonance  with  the  microwave  excitation.  An 
exponential  fit  to  the  data  implies  T\  =  1.3  /is  (line). 

ng  is  then  shifted  to  ng  =  0.171  in  30  ns,  slowly  enough  to  be  adiabatic  but  much  faster  than 
T\ .  The  microwave  excitation  then  no  longer  resonantly  couples  the  ground  and  excited 
state,  and  the  probability  of  being  in  the  excited  state  decays  in  a  time  Tx .  By  averaging 
many  of  the  transient  responses  to  this  stimulus,  we  find  T\  =  1.3  /rs  (Fig.  4).  A  similar 
T\  was  found  in  [4]  for  a  Cooper-pair  box  with  much  smaller  Ec  and  operated  at  ng  =  0.5. 
The  lifetime  is  a  quantity  which  is  insensitive  to  slow  drifts  in  n0f /  and  demonstrates  that 
the  TLS,  which  oscillates  T\  x  cnoi  =  6  x  105  times  before  relaxing  into  its  ground  state,  is 
well  decoupled  from  all  other  sources  of  dissipation. 

We  can  compare  this  long  lifetime  with  the  spontaneous  emission  rate  expected  from  the 
quantum  fluctuations  of  a  generic  electromagnetic  environment.  Calculating  the  spontaneous 
emission  rate  using  Fermi’s  golden  rule  gives 

|T  =  i^c^j  (|)  sin2 (°)sv Ooi  =  <Wi/270  (2) 

where  Sy{y)  =  2hv(R,e(Z0 ))  is  the  voltage  spectral  density  of  the  quantum  fluctuations  of  an 
environment  with  an  impedance  Zq  at  frequency  v  and  sind  =  Ej/Tuoqi  [11].  The  quantity 
Cj  is  the  total  capacitance  of  the  box  to  nearby  metal  traces,  including  intentional  coupling 
to  the  gate  lead  and  other  unintended  capacitive  coupling  (Fig.  1).  We  calculate  Ti  for  a  50  fl 


environment  to  be  between  0.25  and  1  /x s,  extracting  Cj  =  45  ±  15  aF  from  an  electrostatic 
simulation  of  the  chip  layout  [11,  12].  We  do  not  claim  to  have  demonstrated  that  the  lifetime 
is  limited  by  spontaneous  emission;  however,  if  there  are  additional  relaxation  processes, 
either  due  to  the  electrometer  or  fluctuations  of  some  microscopic  degree  of  freedom  in  the 
box,  their  influence  is  at  most  comparable  to  that  of  spontaneous  emission  into  a  typical 
(Z o  ~  50  11)  electromagnetic  environment. 

In  these  experiments,  we  demonstrate  that  a  Cooper-pair  box  is  a  coherent  two-level  sys¬ 
tem  with  a  long  excited-state  lifetime.  With  spectroscopy,  we  determine  the  box’s  Hamilto¬ 
nian  and  estimate  the  rate  of  spontaneous  emission  of  the  box  into  a  typical  environment. 
We  measure  an  excited-state  lifetime  of  box  that  is  remarkable  for  two  reasons.  First,  it 
shows  that  a  quantum-coherent  microelectronic  circuit  can  have  a  T\  that  approaches  the 
limit  set  by  spontaneous  emission  of  a  photon  into  the  electromagnetic  environment.  Second, 
it  is  observed  by  resolving,  on  sub-microsecond  time  scales,  the  decay  of  the  excited-state 
charge  signal  while  the  two-level  system  is  continuously  measured.  Given  the  observed  elec¬ 
trometer  sensitivity  of  4  x  1CT5  e/\/Hz,  the  excited-state  lifetime  is  long  enough  that  a  single 
measurement  can  discriminate  between  the  box  in  its  excited  state  and  the  box  in  its  ground 
state.  In  a  coherent  superposition  of  states  the  box  oscillates  6  x  105  times  before  decaying 
to  the  ground  state,  demonstrating  that  the  circuit  is  a  promising  qubit  implementation  if, 
as  in  [4],  the  sources  of  inhomogeneous  broadening  can  be  overcome. 
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1.  Introduction 

Electrical  engineers  and  physicists  are  naturally  very  interested  in  the  noise 
of  circuits,  amplifiers  and  detectors.  This  noise  has  many  origins,  some  of 
which  are  completely  unavoidable.  For  example,  a  dissipative  element  (a 
resistor)  at  finite  temperature  inevitably  generates  Johnson  noise.  Engi¬ 
neers  long  ago  developed  spectrum  analyzers  to  measure  the  intensity  of 
this  noise.  Roughly  speaking,  these  spectrum  analyzers  consist  of  a  resonant 
circuit  to  select  a  particular  frequency  of  interest,  followed  by  an  amplifier 
and  square  law  detector  (e.g.  a  diode  rectifier)  which  measures  the  mean 
square  amplitude  of  the  signal  at  that  frequency. 

With  the  advent  of  very  high  frequency  electronics  operating  at  low 
temperatures,  we  have  entered  a  new  regime  hw  >  ksT,  where  quantum 
mechanics  plays  an  important  role  and  one  has  to  begin  to  think  about 
quantum  noise  and  quantum-limited  amplifiers  and  detectors.  This  topic  is 
well-studied  in  the  quantum  optics  community  and  is  also  commonplace  in 
the  radio  astronomy  community.  It  has  recently  become  of  importance  in 
connection  with  quantum  computation  and  the  construction  of  mesoscopic 
electrical  circuits  which  act  like  artificial  atoms  with  quantized  energy 
levels.  It  is  also  important  for  understanding  the  quantum  measurement 
process  in  mesoscopic  systems. 

In  a  classical  picture,  the  intensity  of  Johnson  noise  from  a  resistor  van¬ 
ishes  linearly  with  temperature  because  thermal  fluctuations  of  the  charge 
carriers  cease  at  zero  temperature.  One  knows  from  quantum  mechanics, 
however,  that  there  are  quantum  fluctuations  even  at  zero  temperature, 
due  to  zero-point  motion.  Zero-point  motion  is  a  notion  from  quantum  me¬ 
chanics  that  is  frequently  misunderstood.  One  might  wonder,  for  example, 
whether  it  is  physically  possible  to  use  a  spectrum  analyzer  to  detect  the 
zero-point  motion.  The  answer  is  quite  definitely  yes,  if  we  use  a  quantum 
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system!  Consider  for  example  a  hydrogen  atom  in  the  2p  excited  state  lying 
3/4  of  a  Rydberg  above  the  Is  ground  state.  We  know  that  this  state  is 
unstable  and  has  a  lifetime  of  only  about  1  ns  before  it  decays  to  the  ground 
state  and  emits  an  ultraviolet  photon.  This  spontaneous  decay  is  a  natural 
consequence  of  the  zero-point  motion  of  the  electromagnetic  fields  in  the 
vacuum  surrounding  the  atom.  In  fact,  the  rate  of  spontaneous  decay  gives 
a  simple  way  in  which  to  measure  this  zero  point  motion  of  the  vacuum. 
Placing  the  atom  in  a  resonant  cavity  can  modify  the  strength  of  the  noise 
at  the  transition  frequency,  and  this  effect  can  be  measured  via  a  change 
in  the  decay  rate. 

At  finite  temperature,  the  vacuum  will  contain  blackbody  photons  which 
will  increase  the  rate  of  decay  due  to  stimulated  emission  and  also  cause 
transitions  in  the  reverse  direction,  Is  — >  2p,  by  photon  absorption.  With 
these  ideas  in  mind,  it  is  now  possible  to  see  how  to  build  a  quantum 
spectrum  analyzer. 

The  remainder  of  this  article  is  organized  as  follows.  First  we  describe 
the  general  concept  of  a  two-level  system  as  a  quantum  spectrum  analyzer. 
We  next  review  the  Caldeira-Leggett,  formalism  for  the  modelling  of  a 
dissipative  circuit  element,  such  as  a  resistor,  and  its  associated  quantum 
noise.  Then,  a  brief  discussion  of  the  single  Cooper-pair  box,  a  circuit  which 
behaves  as  a  two-level  system  or  qubit,  is  given.  We  then  discuss  the  effects 
of  a  dissipative  electromagnetic  environment  on  the  box,  and  treat  the  case 
of  a  simple  linear,  but  nonequilibrium  environment,  consisting  of  a  classical 
tunnel  junction  which  produces  shot  noise  under  bias.  Finally,  we  describe 
a  theoretical  technique  for  calculating  the  properties  of  a  Cooper-pair  box 
coupled  to  a  measurement  system,  which  will  be  a  nonlinear,  nonequilibrium 
device,  such  as  a  single-electron  transistor.  Equivalently,  this  allows  one 
to  calculate  the  full  quantum  noise  spectrum  of  the  measurement  device. 
Results  of  this  calculation  for  the  case  of  a  normal  SET  are  presented. 


2.  Two- level  systems  as  spectrum  analyzers 


Consider  a  quantum  system  (atom  or  electrical  circuit)  which  has  its  two 
lowest  energy  levels  eo  and  e±  separated  by  energy  Eq\  =  Hloqi.  We  suppose 
for  simplicity  that  all  the  other  levels  are  far  away  in  energy  and  can  be 
ignored.  The  states  of  any  two-level  system  can  be  mapped  onto  the  states 
of  a  fictitious  spin- 1/2  particle  since  such  a  spin  also  has  only  two  states 
in  its  Hilbert  space.  With  spin  up  representing  the  ground  state  (\g))  and 
spin  down  representing  the  excited  state  (|e)),  the  Hamiltonian  is  (taking 
the  zero  of  energy  to  be  the  center  of  gravity  of  the  two  levels) 


H0  =  - 


frw0i 

- <7 


2 


(1) 
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In  keeping  with  the  discussion  above,  our  goal  is  to  see  how  the  rate  of  ‘spin- 
flip’  transitions  induced  by  an  external  noise  source  can  be  used  to  analyze 
the  spectrum  of  that  noise.  Suppose  for  example  that  there  is  a  noise  source 
with  amplitude  f(t)  which  can  cause  transitions  via  the  perturbation1 

V  =  Af  (t)ax,  (2) 

where  A  is  a  coupling  constant.  The  variable  f(t)  represents  the  noise 
source.  We  can  temporarily  pretend  that  /  is  a  classical  variable,  although 
its  quantum  operator  properties  will  be  forced  upon  us  very  soon.  For  now, 
only  our  two-level  spectrum  analyzer  will  be  treated  quantum  mechanically. 

We  assume  that  the  coupling  A  is  under  our  control  and  can  be  made 
small  enough  that  the  noise  can  be  treated  in  lowest  order  perturbation 
theory.  We  take  the  state  of  the  two-level  system  to  be 


In  the  interaction  representation,  first-order  time-dependent  perturbation 
theory  gives 

IM*))  =  IV’(O))  ~^JodT  ^(r)l^(0)).  (4) 

If  we  initially  prepare  the  two-level  system  in  its  ground  state  then  the 
amplitude  to  find  it  in  the  excited  state  at  time  t  is 


iA 

ae  =  ~T. 

[  dr  (e\a x(r)\g) f  (t)  +  0(A2), 

Jo 

(5) 

ii 

i 

1 

[  dr  eiu}°lT f(r)  +  0(A2). 

Jo 

(6) 

We  can  now  compute  the  probability 

pe{t)  =  \ae\2  =  J  j  dTidT2e~l‘JJ01<'Tl~T2)f(Tl)f(T2)  +  0(A3)  (7) 

We  are  actually  only  interested  on  the  average  time  evolution  of  the  system 

Pe(t)  =  dridr2  e--oi(n-r2)  (/(ri)/(T2))  +  0(A 3)  (8) 

1  The  most  general  perturbation  would  also  couple  to  ay  but  we  assume  that  (as  is 
often,  though  not  always,  the  case)  a  spin  coordinate  system  can  be  chosen  so  that  the 
perturbation  only  couples  to  ax ■  Noise  coupled  to  az  commutes  with  the  Hamiltonian 
but  is  nevertheless  important  in  dephasing  coherent  superpositions  of  the  two  states.  We 
will  not  discuss  such  processes  here. 
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We  can  now  perform  a  change  of  variables  in  the  integrals,  r  =  n  —  12  and 
T  =  (n  +  T2)  /2,  and  we  get 

42  ft  fB(T) 

Pe(t)  =  -;  dT  dr  e-^T  (f  (T  +  r/2)/(T  -  r/2))  +  0(JL3)  (9) 

Jo  J —B(T) 

where 

B  (T)  =  T  if  T  <  t/2 

=  t  -  T  if  T  >  t/2. 

Let  us  now  suppose  that  the  noise  correlation  function  is  stationary  (time 
translation  invariant)  and  has  a  finite  but  small  autocorrelation  time  Tf. 
Then  for  t  >  r/  we  can  set  the  bound  B  (T)  to  infinity  in  the  last  integral 
and  write 

A  2  ft  foo 

peW  =  JpJ0  dT  J  dr  e~iw°1T  (/(r)/(0))  +  0(A3)  (10) 

The  integral  over  r  is  effectively  a  sum  of  a  very  large  number  N  ~  t  j  Tf  of 
random  terms  2  and  hence  the  value  undergoes  a  random  walk  as  a  function 
of  time.  Introducing  the  noise  spectral  density 

/+00 

dr  eiuJT  (f  (r)/ (0)),  (11) 

-OO 

we  find  that  the  probability  to  be  in  the  excited  state  increases  linearly 
with  time,3 

Pe(t)  =  Sf(-UJoi)  (12) 

The  time  derivative  of  the  probability  gives  the  transition  rate 

TT  =  woi)  (13) 

Note  that  we  are  taking  in  this  last  expression  the  spectral  density  on  the 
negative  frequency  side.  If  /  were  a  strictly  classical  source  (/(r)/( 0))  would 

2  The  size  of  these  random  terms  depends  on  the  variance  of  /  and  on  the  value  of 
uoiTf  For  uioiTf  ~S>  1  the  size  will  be  strongly  reduced  by  the  rapid  phase  oscillations  of 
the  exponential  in  the  integrand. 

3  Note  that  for  very  long  times,  where  there  is  a  significant  depletion  of  the  probability 
of  being  in  the  initial  state,  first-order  perturbation  theory  becomes  invalid.  However,  for 
sufficiently  small  A,  there  is  a  wide  range  of  times  r/<f  <1/T  for  which  Eq.  12  is  valid. 
Eqs.  13  and  14  then  yield  well-defined  rates  which  can  be  used  in  a  master  equation  to 
describe  the  full  dynamics  including  long  times. 
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be  real  and  Sf(— woi)  =  Sf  (+uoi).  However,  because  as  we  discuss  below 
/  is  actually  an  operator  acting  on  the  environmental  degrees  of  freedom, 
[/(t),/( 0)]  /  0  and  Sf  (-w0i)  /  S/(+w0i)- 

Another  possible  experiment  is  to  prepare  the  two-level  system  in  its 
excited  state  and  look  at  the  rate  of  decay  into  the  ground  state.  The 
algebra  is  identical  to  that  above  except  that  the  sign  of  the  frequency  is 
reversed: 

U  =  jpSf(+w  (14) 

We  now  see  that  our  two-level  system  does  indeed  act  as  a  quantum  spec¬ 
trum  analyzer  for  the  noise.  Operationally,  we  prepare  the  system  either 
in  its  ground  state  or  in  its  excited  state,  weakly  couple  it  to  the  noise 
source,  and  after  an  appropriate  interval  of  time  (satisfying  the  above  in¬ 
equalities)  simply  measure  whether  the  system  is  now  in  its  excited  state  or 
ground  state.  Repeating  this  protocol  over  and  over  again,  we  can  find  the 
probability  of  making  a  transition,  and  thereby  infer  the  rate  and  hence  the 
noise  spectral  density  at  positive  and  negative  frequencies.  Note  that  in  con¬ 
trast  with  a  classical  spectrum  analyzer,  we  can  separate  the  noise  spectral 
density  at  positive  and  negative  frequencies  from  each  other  since  we  can 
separately  measure  the  downward  and  upward  transition  rates.  Negative 
frequency  noise  transfers  energy  from  the  noise  source  to  the  spectrometer. 
That  is,  it  represents  energy  emitted  by  the  noise  source.  Positive  frequency 
noise  transfers  energy  from  the  spectrometer  to  the  noise  source .  In  order  to 
exhibit  frequency  resolution,  A uo,  adequate  to  distinguish  these  two  cases, 
it  is  crucial  that  the  two-level  quantum  spectrometer  have  sufficient  phase 
coherence  so  that  the  linewidth  of  the  transitions  satisfies  the  condition 
loqi/Auj  >  max[kBT/huJoi,  1]. 

In  thermodynamic  equilibrium,  the  transition  rates  must  obey  detailed 
balance  T^/Tf  =  e^^01  in  order  to  give  the  correct  equilibrium  occupancies 
of  the  two  states  of  the  spectrometer.  This  implies  that  the  spectral  densities 
obey 

Sf  (+uoi)  =  e^Sf{-eo  01).  (15) 

Without  the  crucial  distinction  between  positive  and  negative  frequencies, 
and  the  resulting  difference  in  rates,  one  always  finds  that  our  two  level 

4  Unfortunately,  there  are  several  conventions  in  existence  for  describing  the  noise 
spectral  density.  It  is  common  in  engineering  contexts  to  use  the  phrase  ‘spectral  density’ 
to  mean  Sf(+tv)  +  Sf(—u>).  This  is  convenient  in  classical  problems  where  the  two 
are  equal.  In  quantum  contexts,  one  sometimes  sees  the  asymmetric  part  of  the  noise 
S,/(+w)  —  referred  to  as  the  ‘quantum  noise.’  We  feel  it  is  simpler  and  clearer  to 

simply  discuss  the  spectral  density  for  positive  and  negative  frequencies  separately ,  since 
they  have  simple  physical  interpretations  and  directly  relate  to  measurable  quantities. 
This  convention  is  especially  useful  in  non-equilibrium  situations  where  there  is  no  simple 
relation  between  the  spectral  densities  at  positive  and  negative  frequencies. 


rjsdraftl2.tex;  11/10/2002;  0:40;  p.5 


6 


systems  are  completely  unpolarized.  If,  however,  the  noise  source  is  an 
amplifier  or  detector  biased  to  be  out  of  equilibrium,  no  general  relation 
holds. 

We  now  rigorously  treat  the  quantity  /(r)  as  quantum  operator  in  the 
Hilbert  space  of  the  noise  source.  The  previous  derivation  is  unchanged,  and 
Eqs.  (13,14)  are  still  valid  provided  that  we  interpret  the  angular  brackets 
in  Eq.  (8)  as  representing  the  quantum  statistical  expectation  value  for  the 
operator  correlation  (in  the  absence  of  the  coupling  to  the  spectrometer) 

/  +  0O 

dr  eluJT  ^2  Paa  (a\f(,T)\'~f)  <7|/(0)|a)  (16) 

a, 7 

where  for  simplicity  we  have  assumed  that  (in  the  absence  of  the  coupling 
to  the  spectrometer)  the  density  matrix  is  diagonal  in  the  energy  eigen- 
basis  and  time-independent  (but  not  necessarily  given  by  the  equilibrium 
expression).  This  yields  the  standard  quantum  mechanical  expression  for 
the  spectral  density 

/+oo 

a, 7 

=  2irh^2paQ  |(a|/|7)|25(e7  —  eQ  -  hw). 

a, 7 

Substitution  of  this  into  Eqs.  (13,14)  we  derive  the  familiar  Fermi  Golden 
Rule  expressions  for  the  two  transition  rates. 

In  standard  courses,  one  is  not  normally  taught  that  the  transition 
rate  of  a  discrete  state  into  a  continuum  as  described  by  Fermi’s  Golden 
Rule  can  (and  indeed  should!)  be  viewed  as  resulting  from  the  continuum 
acting  as  a  quantum  noise  source.  The  above  derivation  hopefully  provides 
a  motivation  for  this  interpretation. 

One  standard  model  for  the  continuum  is  an  infinite  collection  of  har¬ 
monic  oscillators.  The  electromagnetic  continuum  in  the  hydrogen  atom 
case  mentioned  above  is  a  prototypical  example.  The  vacuum  electric  field 
noise  coupling  to  the  hydrogen  atom  has  an  extremely  short  autocorre¬ 
lation  time  because  the  range  of  mode  frequencies  ua  (over  which  the 
dipole  matrix  element  coupling  the  atom  to  the  mode  electric  field  Ea 
is  significant)  is  extremely  large,  ranging  from  many  times  smaller  than  the 
transition  frequency  to  many  times  larger.  Thus  the  autocorrelation  time  of 
the  vacuum  electric  field  noise  is  considerably  less  than  10~15s,  whereas  the 
decay  time  of  the  hydrogen  2p  state  is  about  lCP9s.  Hence  the  inequalities 
needed  for  the  validity  of  our  expressions  are  very  easily  satisfied. 

Of  course  in  the  final  expression  for  the  transition  rate,  energy  con¬ 
servation  means  that  only  the  spectral  density  at  the  transition  frequency 


(17) 

(18) 
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enters.  However,  in  order  for  the  expression  to  be  valid  (and  in  order  for 
the  transition  rate  to  be  time  independent),  it  is  essential  that  there  be  a 
wide  range  of  available  photon  frequencies  so  that  the  vacuum  noise  has  an 
autocorrelation  time  much  shorter  than  the  inverse  of  the  transition  rate. 

3.  Quantum  Noise  from  a  Resistor 

Instead  of  an  atom  in  free  space,  we  might  consider  a  quantum  bit  capaci- 
tively  coupled  to  a  transmission  line.  The  transmission  line  is  characterized 
by  an  inductance  per  unit  length  i  and  capacitance  per  unit  length  c.  A 
semi-infinite  transmission  line  presents  a  frequency-independent  impedance 
Z  =  Rq  =  \J~tfc  at  its  end  and  hence  acts  like  an  ideal  resistor.  The 
dissipation  is  caused  by  the  fact  that  currents  injected  at  one  end  launch 
waves  which  travel  off  to  infinity  and  do  not  return.  Very  conveniently, 
however,  the  system  is  simply  a  large  collection  of  harmonic  oscillators  (the 
normal  modes)  and  hence  can  be  readily  quantized.  This  representation  of 
a  physical  resistor  is  essentially  the  one  used  by  Caldeira  and  Leggett  [1] 
in  their  seminal  studies  of  the  effects  of  dissipation  on  tunneling.  The  only 
difference  between  this  model  and  the  vacuum  fluctuations  in  free  space 
discussed  above  is  that  the  relativistic  bosons  travel  in  one  dimension  and 
do  not  carry  a  polarization  label.  This  changes  the  density  of  states  as  a 
function  of  frequency,  but  has  no  other  essential  effect. 

The  Lagrangian  for  the  system  is 

r°°  p  1 

C  =  L  iX2j  y/’  (19) 

where  j(x,  t)  is  the  local  current  density  and  q(x,  t )  is  the  local  charge  den¬ 
sity.  Charge  conservation  connects  these  two  quantities  via  the  constraint 

dxj(x,  t )  +  d tq{x,  t )  =  0.  (20) 

We  can  solve  this  constraint  by  defining  a  new  variable 

pX 

0(x,t)  =  /  dx'q(x',t )  (21) 

Jo 

in  terms  of  which  the  current  density  is  j(x,  t )  =  —dtd(x,  t)  and  the  charge 
density  is  q(x,t)  =  dx0(x,t).  For  any  well-behaved  function  9(x,t),  the 
continuity  equation  is  automatically  satisfied  so  there  are  no  dynamical 
constraints  on  the  9  field.  In  terms  of  this  field  the  Lagrangian  becomes 

roo  p  i 

C=  /  dx  -  {dtOf  -  —  (dx0f  (22) 

JO  ^  ^ C 
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The  Euler-Lagrange  equation  for  this  Lagrangian  is  simply  the  wave  equa¬ 
tion  v2d20  —  d20  =  0  where  the  mode  velocity  is  v  =  1  /\/Ic. 

/From  Eq.  (21)  we  can  deduce  that  the  proper  boundary  conditions 
(in  the  absence  of  any  coupling  to  the  qubit)  for  the  6  field  are  9(0,  t)  = 
9(L,t)  =  0.  (We  have  temporarily  made  the  transmission  line  have  a  fi¬ 
nite  length  L .)  The  normal  mode  expansion  that  satisfies  these  boundary 
conditions  is 

9(x,t)  =  \f^^2^Pn(t)  sin  (23) 

where  ipn  is  the  normal  coordinate  and  kn  =  ^ .  Substitution  of  this  form 
into  the  Lagrangian  and  carrying  out  the  spatial  integration  yields  a  set  of 
independent  harmonic  oscillators  representing  the  normal  modes. 


°°  p  1 

£  =  yZn^n)2  ~  -k[ 


^  2 

n=  1 


.2..2 

riPn' 


(24) 


From  this  we  can  find  the  momentum  pn  canonically  conjugate  to  <pn  and 
quantize  the  system  to  obtain  an  expression  for  the  voltage  at  the  end  of  the 
transmission  line  in  terms  of  the  mode  creation  and  destruction  operators 


V  = 'll )  W(».  *)  =  1  (25) 

n=  1 


The  spectral  density  of  voltage  fluctuations  is  then  found  to  be 

2  (“*'1  tQ 

Sv(v)  =  27 rT  -Tr^{n1(hUn)6(u  +  Un)  +  [n1(hQn)  +  l}8(u-fln)},  (26) 

Ju  ‘  ^  Zc 

n=  1 


where  n7(fou)  is  the  Bose  occupancy  factor  for  a  photon  with  energy  frw. 
Taking  the  limit  L  — ►  oo  and  converting  the  summation  to  an  integral  yields 


Sv(  w)  =  2Roh\Lu\{n~/(h\Lu\)@(—Lu)  +  [n^(hw)  +  l]0(u;)},  (27) 

where  0  is  the  step  function.  We  see  immediately  that  at  zero  temperature 
there  is  no  noise  at  negative  frequencies  because  energy  can  not  be  extracted 
from  zero-point  motion.  However  there  remains  noise  at  positive  frequencies 
indicating  that  the  vacuum  is  capable  of  absorbing  energy  from  the  qubit. 

A  more  compact  expression  for  this  ‘two-sided’  spectral  density  of  a 
resistor  is 

2  RqHuj 

^  =  1  _  e-hu>/kBT’ 

which  reduces  to  the  more  familiar  expressions  in  various  limits.  For  ex¬ 
ample,  in  the  classical  limit  /crT  hu>  the  spectral  density  is  equal  to  the 
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Johnson  noise  result5 

Sy(u>)  =  (29) 

which  is  frequency  independent,  and  in  the  quantum  limit  it  reduces  to 

Sy(u)  =  2R0hcoQ(co).  (30) 

Again,  the  step  function  tells  us  that  the  resistor  can  only  absorb  energy, 
not  emit  it,  at  zero  temperature. 

If  we  use  the  engineering  convention  and  add  the  noise  at  positive  and 
negative  frequencies  we  obtain 

Sy(u)  +  Sy(-u>)  =  2-Rofancoth  (31) 

2  kb  1 

for  the  symmetric  part  of  the  noise,  which  appears  in  the  quantum  fluctuation- 
dissipation  theorem[2j.  The  antisymmetric  part  of  the  noise  is  simply 

Sy(u)  -  Sy(~Lo)  =  2RqJtw.  (32) 

This  quantum  treatment  can  also  be  applied  to  any  arbitrary  dissipative 
network[3].  If  we  have  a  more  complex  circuit  containing  capacitors  and 
inductors,  then  in  all  of  the  above  expressions,  Rq  should  be  replaced  by 
R eZ(u>)  where  Z(u)  is  the  complex  impedance  presented  to  the  qubit. 

4.  The  Single  Cooper-Pair  Box:  a  Two-Level  Quantum  Circuit 

The  Cooper-pair  box  (CPB)  is  a  simple  circuit  [4],  consisting  of  a  small 
superconducting  “island” ,  connected  to  a  large  reservoir  via  a  single  small- 
capacitance  Josephson  junction,  depicted  as  a  box  with  a  cross  (Fig.  1).  The 
island  is  charge  biased  by  applying  a  voltage  (Vg)  to  a  nearby  lead,  called 
the  gate,  which  has  a  small  capacitance  to  the  island,  Cg.  The  junction 
is  characterized  by  its  capacitance,  Cj,  and  its  tunnel  resistance,  Rj.  At 
temperatures  well  below  the  transition  temperature  of  the  superconductor 
{Tc  r-j  1.5  K  for  the  usual  Al/AlOx/Al  junctions),  none  of  the  many  (~  109) 
quasiparticle  states  on  the  island  should  be  thermally  occupied,  and  the 
number  of  Cooper-pairs  on  the  island  is  the  only  relevant  degree  of  freedom. 

We  may  then  write  the  Hamiltonian  for  the  box  in  terms  of  the  states  of 
different  numbers  of  pairs  on  the  island,  which  are  eigenstates  of  the  number 
operator,  h|n)  =  n|n).  The  box  Hamiltonian  consists  of  an  electrostatic 
term,  plus  a  Josephson  term  describing  the  coupling  of  the  island  to  the 

5  Note  again  that  in  the  engineering  convention  this  would  be  SV( oj)  =  ARok^T. 


rjsdraftl2.tex;  11/10/2002;  0:40;  p.9 


10 


Figure  1.  a)  Circuit  diagram  of  Cooper-pair  box.  b)  Pseudo-spin  representation  of  the 
energies  of  Cooper-pair  box.  The  density  matrix  for  the  two  pure  eigenstates  lie  along 
the  total  effective  field,  collinear  with  the  z’  axis. 


lead, 


H  —  Heieciro static  “1“  HjosepfiSOn 


AEc^n-ngf\n){n\  -  ^{\n  +  1) {n\  +  h.c.) 


(33) 

(34) 


The  energy  scale  for  the  electrostatic  interaction  is  given  by  the  charging 
energy,  Eq  =  e2/2Cs,  where  Cs  =  C3  +  Cg  is  the  total  island  capacitance, 
while  the  Josephson  energy,  Ej,  is  set  by  the  tunnel  resistance  and  the  gap 
of  the  superconductor, 


hA  A  Rk 
J  8 e2Rj  8  Rj 


(35) 


The  electrostatic  term  is  easily  modulated  by  changing  the  voltage  on  the 
gate;  the  quantity  ng  =  CgVg/2e  that  appears  in  the  Hamiltonian  corre¬ 
sponds  to  the  total  polarization  charge  (in  units  of  Cooper  pairs)  injected 
into  the  island  by  the  voltage  source. 

This  Hamiltonian  leads  to  particularly  simple  behavior  in  the  charge 
regime,  when  the  electrostatic  energy  dominates  over  the  Josephson  cou¬ 
pling,  AEc  3>  Ej.  In  this  case  we  can  restrict  the  discussion  to  only  two 
charge  states,  | n  =  0)  and  |n  =  1).  For  convenience  we  can  reference  the 
energies  of  the  two  states  to  their  midpoint,  Emtj  =  4Ec(l  —  2n9)2,  so  that 
the  Hamiltonian  now  becomes 


-Eei  -Ej  \ 
-Ej  Eei  J 


(36) 


where  Eei  is  the  electrostatic  energy  that  is  now  linear  in  the  gate  charge, 
Eei  =  ^Ec (l—2ng).  It  is  also  now  apparent  that  the  Hamiltonian  is  identical 
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Figure  2.  Energies  (left)  of  ground  and  excited  states  of  a  Cooper-pair  box  with 
Ec  =  Ej  vs.  dimensionless  gate  charge,  ng  =  CgVg/2e.  The  expectation  value  of  a 
charge  measurement,  (n),  (right)  for  the  ground  (solid  line)  and  excited  (dotted  line) 
states  vs.  ng. 


to  that  of  a  fictitious  spin-1/2  particle, 


ZJ  -^el  J  (nr7\ 

H  = - —CTZ  -  —ax,  (37) 

under  the  influence  of  two  psuedo-magnetic  fields,  Bz  =  Eei  and  Bx  =  Ej, 
as  depicted  in  Fig.  1.  In  other  words,  the  box  is  a  qubit  or  two-level  system6. 
The  state  of  the  system  is  in  general  a  linear  combination  of  the  states 
| n  =  0)  and  | n  =  1).  The  state  can  be  depicted  using  the  density  matrix, 
which  corresponds  to  a  point  on  the  Bloch  sphere,  where  the  north  pole 
(+z-direction)  corresponds  to  | n  =  0).  The  ground  and  excited  states  of 
the  system  will  be  aligned  and  anti-aligned  with  the  total  fictitious  field, 
i.e.  in  the  ±z'  directions. 

It  is  also  apparent  from  this  discussion  that  the  states  of  the  box  can 
be  easily  manipulated  by  changing  the  gate  voltage.  The  energies  of  the 
ground  and  excited  states,  as  a  function  of  ng,  are  displayed  in  Figure  2. 
The  energy  difference  between  the  ground  and  excited  bands  varies  from 
AEc  at  ng  =  0, 1,  to  a  minimum  at  the  charge  degeneracy  point,  ng  =  1/2. 
At  this  point,  the  Josephson  coupling  leads  to  an  avoided  crossing,  and  the 
splitting  is  Ej. 

Also  plotted  is  the  expectation  value  of  the  number  operator,  (|h|), 
which  is  proportional  to  the  total  charge  on  the  island.  In  the  geometrical 

6  Of  course,  this  is  an  approximation,  as  there  are  other  charge  states  (|n  =  2),  etc.) 
which  are  possible,  but  require  much  higher  energy.  Even  outside  the  charge  regime  (i.e. 
when  Ej  >  4 Ec)  the  two  lowest  levels  of  the  box  can  be  used  to  realize  a  qubit  [5].  In 
this  case,  the  two  states  do  not  exactly  correspond  with  eigenstates  of  charge,  and  matrix 
elements  are  more  complicated  to  calculate.  Nonetheless,  this  regime  can  also  be  used  as 
an  electrical  quantum  spectrum  analyzer. 
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picture  of  Fig.  lb,  a  measurement  of  charge  (h)  is  equivalent  to  projecting 
the  state  on  the  z-axis,  h  =  ^(1  —  az).  We  see  that  as  the  gate  charge 
is  changed  from  0  to  1,  the  ground  state  is  initially  |n  =  0),  and  the 
character  of  the  ground  and  excited  states  interchange  on  passing  through 
the  degeneracy  point,  leading  to  the  transition  between  (h)  =  0  and  1,  which 
is  broadened  by  quantum  fluctuations  (the  ax  coupling).  At  the  degeneracy 
point,  the  ground  and  excited  states  lie  in  the  ±x-directions,  i.e.  they  are 
symmetric  and  antisymmetric  combinations  of  the  two  charge  states.  In 
general,  we  will  denote  the  ground  and  excited  state  of  the  CPB  at  a 
particular  gate  voltage  as  |  g)  and  |e),  which  are  given  in  terms  of  the  charge 
states  by  \g)  =  cos(0/2)|O)+sin(0/2)|l)  and  |e)  =  —  sin(0/2)|O)+cos(6)/2)|l) 
respectively,  where  0  =  arctan [Ej/Eei]  is  a  function  of  the  gate  voltage. 

A  nice  property  of  the  CPB  in  this  regime  is  that  the  various  matrix 
elements  can  be  calculated  in  a  straightforward  way.  For  example,  the 
expectation  value  of  h  in  the  ground  state,  (g\h\g),  is  therefore  equal  to 
1/2(1  —  (g\az\g))  =  sin2(0/2),  from  which  we  can  find  the  ground  state 
charge  as  shown  in  Fig.  2.  A  perturbation  in  the  gate  charge,  due  for 
example  to  a  fluctuation  or  change  in  the  applied  gate  voltage,  will  lead 
to  a  proportional  change  in  the  electrostatic  energy,  or  the  z-component  of 
the  fictitious  magnetic  held.  Such  a  perturbation  will  cause  both  dephasing 
and  transitions  between  states. 

5.  General  Discussion  of  CPB  Coupled  to  a  Dissipative  Environment 

In  the  previous  section  we  described  the  Hamiltonian  and  the  eigenstates  for 
a  Cooper-pair  box  which  is  “charge-biased,”  i.e.  controlled  with  a  voltage 
applied  to  a  gate  capacitor  Cg,  as  shown  in  Fig.  1.  In  our  earlier  treatment 
of  the  box,  the  voltage  and  the  dimensionless  gate  charge,  ng  were  treated 
as  fixed  parameters  of  the  Hamiltonian  (c-numbers).  In  this  case,  the  box’s 
evolution  is  purely  deterministic  and  conservative.  However,  it  is  impossible, 
even  in  principle,  to  control  such  a  voltage  with  arbitrary  precision  at  all 
frequencies.  In  Fig.  1,  the  idealized  source  of  the  gate  voltage  is  drawn  in 
series  with  an  impedance  Z(u))  of  the  gate  lead.  Generally,  this  gate  lead 
will  be  connected  to  external  wiring  (a  transmission  line),  with  a  typical 
real  impedance  comparable  to  the  impedance  of  free  space  (~  50  H)  at 
the  microwave  transition  frequencies  of  the  box.  /From  the  fluctuation- 
dissipation  theorem  we  know  that  this  impedance  will  introduce  noise  on 
the  gate  voltage,  even  at  zero  temperature. 

There  are  several  effects  of  the  voltage  noise  on  the  box,  or  the  coupling 
of  our  spin-1/2  circuit  to  the  many  external  degrees  of  freedom  represented 
by  the  gate  impedance.  First,  even  at  zero  temperature,  we  will  find  a  finite 
excited  state  lifetime,  Tl,  for  the  box.  Second,  at  finite  temperature,  we  will 
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find  a  finite  polarization  of  our  psuedo-spin,  i.e.  some  steady-state  proba¬ 
bility  to  find  the  spin  in  its  excited  state.  Finally,  the  gate  noise  introduces 
a  random  effective  field  felt  by  the  spin,  and  a  loss  of  phase  coherence  for  a 
superpostition  state.  It  is  this  last  effect  which  is  most  important  in  making 
high-fidelity  qubits  and  performing  quantum  computations,  but  it  is  the 
first  two  which  depend  most  explicitly  on  the  quantum  nature  of  the  noise. 
We  deal  in  this  manuscript  with  only  these  first  two  features  of  the  box’s 
coupling  to  the  electromagnetic  environment,  and  ignore  the  dephasing7.  Of 
course,  the  other  parameter  in  the  Hamiltonian,  the  Josephson  energy,  can 
in  principle  fluctuate,  especially  as  in  many  experiments  the  box’s  junction 
is  split  into  a  small  SQUID  in  order  to  provide  external  tuning  of  Ej  with 
an  applied  flux.  We  concentrate  here  only  on  the  voltage  noise  (fluctuations 
in  the  az  part  of  the  Hamiltonian)  for  simplicity. 

We  begin  with  a  very  simple  treatment  of  the  dynamics  of  the  two-level 
system  under  the  influence  of  the  gate  voltage  noise.  We  are  interested  in 
the  ensemble-averaged  behavior  of  our  psuedo-spin,  which  is  best  described 
using  the  density  matrix  approach,  and  is  detailed  in  Section  7  on  the 
coupling  of  the  box  to  a  measuring  SET.  The  basic  effects,  apart  from 
dephasing,  however,  can  be  captured  simply  by  examining  the  probabilities 
pg  and  pe  of  finding  the  box  in  its  ground  (|<?))  or  excited  (|e))  states. 
The  noise  of  the  external  environment  can  drive  transitions  from  ground  to 
excited  state  and  vice-versa,  at  rates  Tj  and  Tj  ,  respectively.  The  coupled 
master  equations  for  these  probabilities  are 


%  =P9T]~Pe^l 


(38) 

(39) 


Of  course  conservation  of  probability  tells  us  that  pe  +  pg  =  1,  so  we 
introduce  the  polarization  of  the  spin-1/2  system,  P  =  pg  —  pe.  In  steady- 
state,  the  detailed  balance  condition  is  peT j_  =  pgY p  The  two  rates  Tj 
and  Ti  are  related  by  Equations  13  and  14  to  the  spectral  densities  of  the 
noise  at  negative  and  positive  frequencies.  We  see  immediately  that  if  the 
spectral  density  is  symmetric  (classical!),  then  the  rates  for  transitions  up 
and  down  are  equal,  the  occupancies  of  the  two  states  are  exactly  equal,  and 
the  polarization  of  the  psuedo-spin  is  identically  zero.  It  is  the  quantum,  or 
antisymmetric,  part  of  the  noise  which  gives  the  finite  polarization  of  the 
spin.  Even  in  NMR,  where  the  temperature  is  large  compared  to  the  level 
splitting  {JvjJqi  <  ksT),  this  effect  is  well-known  and  crucial,  as  the  small 
but  non-zero  polarization  is  the  subject  of  the  field! 


'  For  a  nice  recent  treatment  of  dephasing  in  .Josephson  junctions,  see  Ref.  [6] . 
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Solving  for  the  steady-state  polarization,  we  find 

p  =  rl  ~  rT  =  ^(+^01 )  ~  S{-um)  ,  > 

ss  ri  +  rT  s{+uQ1)  +  s(-uoi)  y  } 

An  measurement  of  the  steady-state  polarization  allows  one  to  observe  the 
amount  of  asymmetry  in  the  noise,  so  the  two-level  system  is  a  quantum 
spectrum  analyzer. 

If  we  can  create  a  non-equilibrium  polarization,  P  =  Pss  +  A P  (a  pure 
state  is  not  necessary)  of  our  two-level  system,  we  expect  it  to  return  to  the 
steady  state  value.  Substituting  the  modified  probabilities  pe(t)  =  pess  — 
AP(t)/2  and  pg(t )  =  p9ss  +  AP(t)/2  into  our  master  equations  above,  we 
find  an  equation  for  the  deviation  of  the  polarization 


d(AP(t)) 

dt 


-AP^^  +  rj. 


(41) 


Thus  the  system  decays  to  its  steady-state  polarization  with  the  relaxation 
rate  Ti  =  Tj  +  Tj  =  (A/h)2[S(—co oi)  +  <S'(+woi)]  related  to  the  total 
noise  at  both  positive  and  negative  frequencies.  In  NMR,  the  time  1/Ti 
is  referred  to  as  Tl.  In  the  zero-temperature  limit,  there  is  no  possibility  of 
the  qubit  absorbing  energy  from  the  environment,  so  Tj  =  0,  and  we  find 
full  polarization  P  =  1,  and  a  decay  of  any  excited  state  population  at  a 
rate  Tj  =  1/TI  which  is  the  spontaneous  emission  rate. 

It  is  worth  emphasizing  that  a  quantum  noise  source  is  always  char¬ 
acterized  by  two  numbers  (at  any  frequency),  related  to  the  positive  and 
negative  frequency  spectral  densities,  or  to  the  symmetric  and  antisym¬ 
metric  parts  of  the  noise.  These  two  quantities  have  different  effects  on  a 
two-level  system,  introducing  a  finite  polarization  and  finite  excited-state 
lifetime.  Consequently,  a  measurement  of  both  the  polarization  and  Tl  of  a 
two-level  system  is  needed  to  fully  characterize  the  quantum  noise  coupled 
to  the  qubit.  Such  measurements  in  electrical  systems  are  now  possible, 
and  some  of  us  [7]  have  recently  performed  such  a  characterization  for 
the  specific  case  of  a  CPB  coupled  to  a  superconducting  single-electron 
transistor. 

Our  discussion  in  this  section  uses  the  language  of  NMR  to  describe 
the  effects  on  the  two-level  system.  There  are,  however,  several  possible 
protocols8  for  measuring  the  quantum  noise,  and  several  different  “basis 
sets”  or  measured  quantities  which  describe  the  noise  process  or  the  quan¬ 
tum  reservoir  to  which  the  two-level  system  is  coupled.  Table  1  contains 
a  “translation”  between  the  specific  pairs  of  quantities  that  are  commonly 


8  The  idea  of  watching  the  decay  from  the  pure  states  |e)  and  | g)  to  measure 
SV(±|a>oi|)  separately  was  described  in  Section  2. 
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TABLE  I.  Different  ways  to  characterize  a  quantum  reservoir. 


Fermi  Golden  Rule 

rr(w)  —  jtSv(+ M) 

FjM  =  j?Sv{-\w\) 

Fluct.-Diss.  Relation 

n-y(u>)t=  2rT/(r|  —  rT) 

Rc[^(w)]  =  A(ri-rT) 

NMR 

7T  ^  (r,  -  F-)  1 

p  =  (ri-rT)/(ri  +  rT) 

Quantum  Optics 

B  Einstein  —  T  J 

A-Einstein  = 

used  in  different  disciplines,  and  their  relation  to  the  positive  and  negative 
noise  spectral  densities.  In  all  cases,  though,  two  separate  numbers  are 
required  to  specify  the  properties  of  a  quantum  reservoir. 


6.  The  Box  Coupled  to  an  Ohmic  Environment 


We  can  now  proceed  to  the  effects  of  a  specific  dissipative  coupling  to  the 
Cooper-pair  box.  The  noise  on  the  gate  voltage  will  lead  to  a  fluctuation  of 
the  gate  charge  parameter,  ng,  and  thus  to  a  fluctuation  of  the  electrostatic 
energy,  i.e.  the  az  term  in  the  Hamiltonian  (Eq.  37).  Depending  on  the 
average  value  of  ng,  this  fluctuation  will  consist  of  fluctuations  which  are 
both  longitudinal  (||  to  a'z)  and  transverse  (_L  to  a'z).  To  calculate  the  rates 
of  transitions  between  the  states  |e)  and  \g),  we  need  to  find  the  coupling 
strength  A  of  this  perturbation  in  the  ax  direction.  Referring  to  Fig.  1, 
we  see  that  az  =  cos(0)(t(,  —  sin(0)o-^,.  If  we  let  the  gate  charge  now  be 
ng{t)  =  ng  +  5ng{t ),  we  may  rewrite  the  Hamiltonian  Eq.  37  in  the  new 
eigenbasis  as 

H  =  ~^r(J,z  +  cos(0)5ng(t)az  —  4Ec  sin (9)8ng{t)a'x  .  (42) 

Zjfb 

The  time-varying  term  in  the  a'z  direction  will  effectively  modulate  the 
transition  frequency,  otoi  =  Eqi/H,  and  cause  dephasing.  In  terms  of  the 
gate  voltage  noise,  V(t),  the  a'x  perturbation  term  has  the  form  AV(t)ax  = 
ensva.(9)V {t)a'x,  where  e  is  the  electron’s  charge  and  n  =  Cg/C s  is  the 
capacitive  coupling.  Using  Eq.  14,  we  find 


Fj.  =  (J^)2  k2  sin2 (6) Sv(+u>01).  (43) 

If  the  environment  is  effectively  at  zero  temperature  (tua oi  3>  k^T),  then 
5y(+a;oi)  =  2hwoiR,  and  the  quality  factor  of  the  transition  is 


Q  —  woi  /Tj_ 


1  Rk 

k2  sin2  9  4irR 


(44) 
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where  Rk  =  h/e 2  is  the  resistance  quantum. 

For  a  finite  temperature,  we  have  rates  in  both  directions,  and  the 
polarization  of  the  psuedo-spin  is  given  by  the  ratio  of  the  antisymmetric 
(Eq.  32)  to  symmetric  (Eq.  31)  spectral  densities,  P  =  t&nh.(hjJoi/2k-p,T), 
as  one  expects  for  any  two-level  system  at  temperature  T.  An  example  of 
the  average  charge  state  of  a  Cooper-pair  box  at  finite  temperature,  and  of 
the  polarization  and  equilibration  time  Tl,  are  shown  in  Figure  3.  As  the 
gate  voltage  is  varied,  the  transition  frequency  of  the  box  changes  from  a 
maximum  near  ng  =  0,  to  a  minimum  ujq\  =  Ej/h  at  the  degeneracy  point 
ng  =  0.5  and  then  back  again.  We  see  that  the  states  of  the  box  are  generally 
most  “fragile”  near  the  avoided  crossing.  First,  the  energy  splitting  is  a  min¬ 
imum  here,  leading  to  the  lowest  polarization  of  the  psuedo-spin.  Second, 
because  the  eigenstates  \g)  and  |e)  point  in  the  ax  direction,  the  matrix 
elements  for  the  voltage  fluctuations  of  the  environment  are  maximal,  i.e. 
the  noise  is  orthogonal  to  the  spin.  This  also  implies  that  the  dephasing 
effects  are  minimal  at  this  degeneracy  point,  which  offers  great  advantages 
for  improving  the  decoherence  times  [5],  but  the  excited  state  lifetimes  are 
smallest  at  this  point.  One  also  sees  that  the  lifetimes  become  large  away 
from  the  degeneracy,  where  the  electrostatic  energy  dominates  over  the 
Josephson  energy,  which  offers  advantages  when  measuring  the  charge  state. 
In  the  limit  that  Ej  could  be  suppressed  to  zero  during  a  measurement,  the 
matrix  elements  (for  voltage  noise)  vanish,  and  a  quantum  non-demolition 
(QND)  measurement  [8]  could  be  performed.  The  idea  of  using  the  qubit 
as  a  quantum  spectrum  analyzer  is  precisely  the  reverse,  where  we  measure 
the  “destruction”  in  the  two-level  system  (i.e.  inelastic  transitions  caused 
by  the  coupling  of  the  states  to  the  environment),  in  order  to  learn  about 
the  quantum  noise  spectrum  of  the  fluctuations. 

The  Cooper-pair  box  can  of  course  also  be  used  to  measure  the  more  in¬ 
teresting  spectral  densities  of  nonequilibrium  devices.  The  simplest  example 
is  to  replace  the  gate  resistance  by  a  tunnel  junction.  If  we  arrange  to  bias 
the  junction  using,  e.g.  an  inductor,  a  dc  current  I  and  an  average  dc  voltage 
V  can  be  maintained  across  the  junction.  Classically,  the  current  noise  of 
such  a  tunneling  process  is  frequency  independent,  Sj  =  2 el.  The  voltage 
noise  density  presented  to  the  CPB’s  gate  would  then  be  Sy  =  2 elR^, 
where  Rr  is  the  junction  tunnel  resistance.  In  fact,  this  “white”  spec¬ 
tral  density  can  only  extend  up  to  frequencies  of  order  u  =  eV/h,  the 
maximum  energy  of  electrons  tunneling  through  the  junction.  The  correct 
high-frequency  form  of  the  symmetrized  noise  density  was  calculated  by 
Rogovin  and  Scalapino  [9], 


Sy  (w) 


R  (Hco  +  eV)  coth 


TvjJ  +  eV 
2  ksT 


+  R  (Hu 


eV)  coth 


Huj  —  eV 


2 kBT  J  ’ 
(45) 
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Figure  3.  The  box  coupled  to  an  equilibrium,  Ohmic  environment,  i.e.  a  resistor,  a) 
Two-sided  noise  spectral  density,  of  the  voltage,  Sv{ w),  for  a  resistor  at  T=0  (solid 
line)  and  finite  temperature  (dashed  line)  b)  Average  charge  of  box  with  Ec  =  1A', 
Ej  =  0.5A'  when  coupled  with  strength  n  =  0.01  to  a  resistor  with  resistance  R  =  50  Q 
and  temperature  T  =  0.5K.  c)  Polarization  (dotted  line)  and  relaxation  time  T1  for  the 
same  parameters.  Full  line  is  the  rate  of  spontaneous  emission,  i.e.  T1  at  zero  temperature. 


and  was  indirectly  measured  in  a  mesoscopic  conductor  using  a  conventional 
spectrum  analyzer  [10].  This  noise  can  also  be  expressed  [11]  in  its  two-sided 
form 


Sy  (uj) 


[hjjj  +  eV)  Rt  (hio  —  eV )  Rt 

_ huj  +  eV  _ hu  —  eV  "> 

1  —  e  kBT  1  —  e  ksT 


(46) 


and  is  displayed  in  Fig.  4.  Notice  that  the  antisymmetric  part  of  this  noise  is 
the  same  as  that  of  the  ordinary  resistor,  SV(+u;)  —  <SV(— lu)  =  2 HloRt,  and 
is  independent  of  the  voltage.  Also  shown  in  Fig.  5  is  the  polarization  and 
relaxation  time,  Tl,  of  a  CPB  coupled  to  a  shot  noise  environment.  We  see 
that  full  polarization  is  achieved  only  when  hw oi  3>  eV.  For  low  transition 
frequencies,  the  polarization  is  inversely  proportional  to  the  current  through 
the  junction.  Aguado  and  Kouwenhoven  [11]  have  described  the  use  of  a 
double  quantum  dot  as  a  two-level  system  to  probe  this  behavior  of  the 
shot  noise. 

Given  our  discussion  so  far,  it  is  now  interesting  to  ask  what  the  effects 
of  a  real  quantum  measurement  on  a  quantum  circuit  will  be.  A  quantum 
measurement  device  will  in  general  be  neither  linear,  Ohmic,  nor  equilib- 
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Figure  4-  The  box  coupled  to  an  nonequilibrium,  Ohmic  environment,  i.e.  a  50  SI 
tunnel  junction,  a)  Two-sided  noise  spectral  density,  of  the  voltage,  Sv( w),  for  a  junction 
at  T=0.02  K,  with  zero  voltage  (dotted  line),  and  increasing  bias  voltages  (solid  and 
dashed  lines),  b)  Average  charge  of  box  with  Ec  =  IK,  Ej  =  0.5K  when  coupled  with 
strength  k  =  0.01  to  a  junction  biased  at  eV  =  1.5  K.  c)  Polarization  (dotted  line)  and 
relaxation  time  T1  for  the  same  parameters  (solid)  and  for  T=0,  V=0  (dashed  line). 


rium.  Obviously,  if  we  hope  to  characterize  this  measurement  process,  and 
to  understand  what  one  will  observed  when  the  qubit  is  coupled  to  the 
noise  processes  of  the  measuring  device,  we  will  need  to  calculate  the  full 
quantum  (two-sided!)  noise  spectrum  of  the  amplifier  or  detector. 

7.  Single-Electron  Transistor  Coupled  to  a  Two-Level  System 

We  have  seen  in  previous  sections  that  a  two-level  system  (TLS)  may  be 
used  as  a  “spectrum  analyzer”  to  measure  quantum  noise.  Here,  we  use  this 
technique  to  theoretically  calculate  quantum  noise.  Instead  of  simply  study¬ 
ing  the  “noisy”  system  of  interest  in  isolation,  one  can  study  a  composite 
system  consisting  of  the  “noisy”  system  coupled  to  a  TLS;  by  calculating 
the  relaxation  and  excitation  rates  of  the  TLS,  one  can  efficiently  calcu¬ 
late  the  quantum  noise  of  interest9.  We  demonstrate  the  usefulness  of  this 

9  Note  that  the  spirit  of  our  approach  is  similar  to  that  employed  in  the  theory  of  full 
counting  statistics  [12].  There  too  one  attaches  an  auxiliary  spin  1/2  to  the  scattering 
system  of  interest,  and  studies  the  dynamics  of  the  fully  coupled  system  to  obtain  the 
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technique  by  outlining  a  calculation  of  the  quantum  charge  noise  of  a  single 
electron  transistor  (SET).  This  is  an  important  example,  as  when  an  SET  is 
used  as  an  electrometer,  it  is  this  noise  which  determines  the  measurement 
backaction. 

The  SET  consists  of  a  metallic  island  attached  via  tunnel  junctions  to 
source  and  drain  reservoirs.  It  is  described  by  the  Hamiltonian: 

HsET  =  Y  (£fc  _  clacka  +  Ec(n  ~  N)2  +  HT 
k,a=L,R,I 

Hr  =  t  Y  ( F*CkICqa  +  h-C •) 

k,q,a=L,R 

The  first  term  in  Hset  describes  the  kinetic  energy  of  electrons  in  the  leads 
(' a  =  L,R)  and  on  the  island  ( a  =  I).  The  second  term  is  the  Coulomb 
charging  energy  which  depends  on  n,  the  number  of  excess  electrons  on 
the  island.  This  interaction  term  can  be  tuned  by  changing  the  voltage 
on  a  nearby  gate  electrode  which  is  capacitively  coupled  to  the  island;  M 
represents  the  dimensionless  value  of  this  voltage.  Finally,  Ht  describes  the 
tunneling  of  electrons  through  the  two  SET  tunnel  junctions;  the  conduc¬ 
tance  of  each  junction  (in  units  of  e2/h)  is  given  by  g  =  4n2t2i/Q,  with  i/q 
being  the  density  of  states.  F'  is  an  auxiliary  operator  which  increases  n 
by  one:  [n,  =  F'.  For  simplicity,  we  assume  that  the  two  junctions  of 

the  SET  are  completely  symmetric  (i.e.  equal  junction  capacitances  and 
dimensionless  conductances) . 

Throughout  this  section,  we  will  be  interested  in  the  regime  of  sequential 
tunneling  in  the  SET,  where  transport  involves  energy-conserving  transi¬ 
tions  between  two  charge  states  of  the  SET  island,  say  n  =  0  and  n  =  1. 
These  transitions  are  described  by  simple  rates,  which  can  be  derived  via 
Fermi’s  Golden  rule: 


(47) 

(48) 


■pQ! 

nd=l,n 

7(A  E) 
A  K±l,n 


7([A  E]an±hn) 
gAE/h 

l  _  e—AE/(kgT) 

—  da,n)  fiVns 


(n±  -  -aA  ±  (\ 


V 


(49) 

(50) 

(51) 


T“±i  n  is  the  tunneling  rate  from  the  charge  state  n  to  n  ±  1  through 
junction  cc;  A E  is  the  energy  gained  in  making  the  tunneling  transition, 
and  includes  contributions  both  from  the  drain-source  voltage  Vos  and 
from  the  charging  energy.  Sequential  tunneling  is  the  dominant  transport 
mechanism  when  the  junction  conductances  are  small  (i.e.  g/( 2tt)  <C  1), 

statistics  of  charge  transfer  in  the  scatterer. 


rjsdraftl2.tex;  11/10/2002;  0:40;  p.19 


20 


and  the  dimensionless  gate  voltage  M  is  not  too  far  away  from  a  charge 
degeneracy  point.  Sequential  tunneling  is  the  most  important  regime  for 
measurement  applications,  as  it  yields  the  largest  SET  currents. 

At  low  temperatures,  only  tunnel  events  which  follow  the  voltage  are 
possible.  There  are  thus  there  only  two  relevant  rates:  n  =  0  — »  1  transitions 
occur  through  the  left  junction  at  a  rate  T[q,  while  n  =  1  — ►  0  occur  through 
the  right  junction  at  a  rate  T q^.  The  average  current  will  be  given  by: 


(I)  =eT  =  e 


pL  pR 
1  IQ1  01 

rfe  +  rft 


(52) 


We  are  interested  in  calculating  Sq(uj),  the  quantum  noise  associated 
with  fluctuations  of  the  charge  on  the  central  island  of  the  SET.  It  is  defined 
as: 


dt(n(t)n( 0))e  'lujt 


(53) 


Note  that  we  can  equivalently  think  of  Sq  as  describing  the  voltage  noise 
of  the  island,  as  ViSiand  =  en/Cz,  where  CG  is  the  total  capacitance  of  the 
island.  In  two  limits,  the  form  of  the  island  charge  noise  can  be  anticipated. 
For  uj  — >  0,  the  noise  will  correspond  to  classical  telegraph  noise-  the  island 
charge  n  fluctuates  between  the  values  0  and  1,  with  Poisson-distributed 
waiting  times  determined  by  the  rates  T(q  and  T^.  We  thus  expect  a 
symmetric,  Lorentzian  form  [13]  for  the  noise  at  low  frequencies: 


9F 

s^>^  +  (rj.0  +  ry  (“«^>  <54) 

For  large  frequencies  |o;|  3>  Ec,  we  expect  that  correlations  due  to  the 
charging  energy  will  have  no  influence  on  the  noise.  The  system  will  effec¬ 
tively  look  like  two  tunnel  junctions  in  parallel,  and  we  can  use  the  results 
of  Sec.  6  for  the  corresponding  voltage  noise.  Noting  that  each  junction 
effectively  consists  of  a  resistor  and  capacitor  in  parallel,  we  have  at  zero 
temperature: 


<SqM 


CJ 2  (j2 

=  x  Sv(v)  ->  [2fouRe  Ztot&(u)\ 


UJ 


.<?  4ecv 

2n  h  I 


-0M 


+  UJ2 


(M  >  Ec) 

(55) 


Note  that  Sq(lu)  decays  as  1  /uj  at  large  frequencies,  whereas  Eq.  (54)  for 
classical  telegraph  noise  decays  as  1/a j2. 

Given  these  two  limiting  forms,  the  question  now  becomes  one  of  how 
the  SET  interpolates  between  them.  One  might  expect  that  the  two  results 


rjsdraftl2.tex;  11/10/2002;  0:40;  p.20 


21 


should  simply  be  added  in  quadrature,  but  as  with  combining  thermal  and 
quantum  noise  (see  Sec.  6),  this  is  approach  is  too  simple.  A  completely 
quantum  mechanical  way  of  calculating  the  noise  for  any  frequency  is 
needed.  This  was  recently  accomplished  by  Johansson  et.  al  [14]  using  an 
extension  of  a  technique  developed  by  Scholler  and  Schon  [15].  Here,  we  re¬ 
derive  their  results  using  the  coupled  system  approach  outlined  above.  This 
method  is  physically  motivated  and  allows  for  a  heuristic  interpretation  of 
the  final  result. 

8.  SET  Coupled  to  a  Qubit 

We  now  consider  a  system  where  the  SET  is  coupled  to  a  two-level  system 
(i.e.  a  qubit),  with  a  coupling  Hamiltonian  which  can  induce  transitions 
in  the  TLS.  Using  spin  operators  to  describe  the  qubit,  and  assuming 
operation  at  the  degeneracy  point  for  simplicity,  where  the  transitions  are 
fastest10,  we  have: 

H  =  HSet  -  +  Aazn,  (56) 

where  ft  is  the  qubit  splitting  frequency11,  and  A  is  the  coupling  strength. 
We  can  define  the  rates  T-f  and  Tj  which  are,  respectively,  the  rate  at  which 
the  qubit  is  excited  by  the  SET,  and  the  rate  at  which  the  qubit  is  relaxed 
by  the  SET.  For  a  weak  coupling  ( A  —>  0),  one  has  (c.f.  Eq.  14,13  in  Sec.  2): 

rm  =  ^-sQ(±ft)  (57) 

Eq.  (57)  tells  us  that  if  we  know  the  rates  Tt  and  T  |  for  a  weakly  coupled 
system  at  an  arbitrary  splitting  frequency  ft,  we  know  the  quantum  noise 
SQ(ft)  at  all  frequencies.  This  is  the  essence  of  the  technique  previously 
described,  in  which  a  qubit  acts  as  a  quantum  spectrum  analyzer  of  noise. 
Here,  we  mimic  this  approach  theoretically  by  obtaining  Tt  and  Tj  from  a 
direct  analysis  of  the  coupled  system  in  the  limit  of  weak  coupling  (A  — >  0). 
The  object  of  interest  is  the  reduced  density  matrix  p  which  describes  both 
the  charge  n  of  the  transistor  island  and  the  state  of  the  qubit.  We  are 
interested  in  two  quantities.  First,  what  is  the  stationary  state  of  the  qubit? 
The  stationary  populations  of  the  two  qubit  states  (which  are  determined 

10  This  amounts  to  maximizing  the  “destruction”  due  to  the  SET’s  noise,  and  the 
case  where  9  =  tt/2,  the  qubit  eigenstates  are  in  the  az  =  ax  direction,  and  the  SET’s 
perturbation  is  in  the  —a'x  =  az  direction  (c.f.  Eq.  42).  After  the  noise  of  the  SET  is 
found,  we  can  then  recalculate  the  effects  on  the  qubit  at  various  ng  or  values  of  9  by 
including  the  modified  matrix  elements  in  the  coupling  coefficient,  A. 

11  Henceforth  we  use  Q.  for  the  transition  frequency,  instead  of  the  previous  notation 
woi  =  Eoi/h,  for  compactness. 
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from  the  time-independent  solution  for  p)  will  tell  us  the  polarization  of 
the  qubit,  and  the  amount  of  asymmetry  in  the  noise  (c.f.  Eq.40).  Second, 
how  quickly  do  the  qubit  populations  relax  to  their  stationary  value?  This 
relaxation  will  be  described  by  a  time-dependent  solution  of  p  characterized 
by  a  mixing  rate  TmiX  which  is  the  sum  of  Tj  and  Tj  (c.f.  Eq.  41).  From 
these  two  results  we  can  solve  for  the  individual  values  of  T-|y ^ . 

To  deal  with  the  dynamics  of  p,  we  make  use  of  the  fact  that  sequential 
tunneling  processes  are  completely  described  by  lowest-order  perturbation 
theory  in  the  tunneling  Hamiltonian  Hp .  Keeping  only  second  order  terms 
(there  are  no  non- vanishing  first  order  terms),  one  obtains  the  following 
standard  evolution  equation  in  the  interaction  picture: 

^tP^  =  ~~hJ  dt,(XHT(t),  pF]])  (58) 

The  angular  brackets  denote  the  trace  over  the  single-particle  degrees  of 
freedom  in  the  SET  leads  and  island;  pp  is  the  equilibrium  density  matrix 
corresponding  to  the  state  of  these  degrees  of  freedom  in  the  absence  of 
tunneling.12  Note  that  a  similar  density  matrix  analysis  of  a  qubit  coupled 
to  a  SET  was  recently  discussed  by  Makhlin  et.  al  [16];  unlike  the  present 
case,  these  authors  restricted  attention  to  a  vanishingly  small  splitting 
frequency  H. 

To  make  progress  with  Eq.  (58),  we  make  a  Markov  approximation, 
which  involves  replacing  pit')  on  the  right-hand  side  with  p(t).  This  is 
permissible  as  we  are  interested  in  the  slow  dynamics  of  p.  We  want  to  find 
both  the  stationary  solution  of  p,  for  which  the  Markov  approximation  is 
exact,  and  the  mixing  mode,  a  mode  whose  time  dependence  is  oc  e-^  T+ri)*. 
This  mode  is  also  arbitrarily  slow  in  the  weak  coupling  limit  A  — ►  0  of 
interest.  Finding  the  stationary  mode  and  the  mixing  mode  correspond  to 
evaluating  the  polarization  and  T1  of  the  qubit  (c.f.  Eq.  40  and  Eq.  41),  as 
was  shown  earlier  for  the  master  equation  of  the  probabilities  in  Section  5. 
Note  that  the  Markov  approximation  should  be  made  in  the  Schrodinger 
picture,  as  it  is  in  the  Schrodinger  picture  that  p  will  be  nearly  stationary 
(i.e.  all  oscillations  associated  with  the  qubit  splitting  frequency  H  will  be 
damped  out  in  the  long-time  limit). 

Evaluation  of  Eq.  (58)  in  the  Markov  approximation  results  in  the  ap¬ 
pearance  of  rates  which  are  generalizations  of  those  given  in  Eq.  (49).  Now, 
however,  these  rates  depend  on  the  initial  and  final  state  of  the  qubit¬ 
tunneling  transitions  can  simultaneously  change  both  the  charge  state  of 
the  SET  island  and  the  state  of  the  qubit.  The  resulting  equation  is  most 

12  In  the  diagrammatic  language  of  Ref.  [15],  Eq.  (58)  is  equivalent  to  keeping  all  ( Ht )2 
terms  in  the  self-energy  of  the  Keldysh  propagator  governing  the  evolution  of  p. 
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easily  presented  if  we  write  the  reduced  density  matrix  p  in  the  basis  of 
eigenstates  at  zero  tunneling.  For  each  value  of  island  charge  n,  there  is  a 
different  qubit  Hamiltonian,  and  correspondingly  a  different  a  qubit  ground 
state  | gn)  and  excited  state  \en).  When  a  tunneling  event  occurs  in  the  SET, 
there  is  a  sudden  change  in  the  qubit  Hamiltonian.  As  the  qubit  ground 
and  excited  states  at  different  values  of  n  are  not  orthogonal,  tunneling 
transitions  in  the  SET  are  able  to  cause  “shake-up”  transitions  in  the  qubit. 
In  the  limit  A  — >  0,  the  relevant  matrix  overlaps  are  given  by: 

...  If  Aim  —  n)  \ 2  ,  .  ,  ,  . 

( 9m\gn )  =  1  —  y  Q  J  =  (em|en)  (59) 

,  .  .  Aim  —  n)  .  . 

(em|  gn)  =  - ^ -  (60) 

Dehning  the  frequency  dependent  rate  rn±i)Tl(u;)  as: 

r'n±l,n(^)  =  £  +  M.  (6D 

a=L,R 

where  A E  and  j(A E)  are  defined  in  Eqs.  (51)  and  (50),  the  required  tunnel 
rates  take  the  form: 

r  m,n  =  rm,„(0)  =  rm,n(±n)  (62) 

The  T+  rates  correspond  to  tunneling  events  where  the  qubit  is  simul¬ 
taneously  relaxed,  and  thus  there  is  an  additional  energy  H  available  for 
tunneling.  For  large  H,  tunneling  processes  which  are  normally  energet¬ 
ically  forbidden  can  occur  if  they  are  accompanied  by  qubit  relaxation. 
Similarly,  the  rates  describe  tunnelling  events  where  the  qubit  is  simul¬ 
taneously  excited,  with  the  consequence  that  there  is  less  energy  available 
for  tunneling. 

Returning  to  the  evolution  equation  Eq.  (58),  note  that  we  do  not 
need  to  track  elements  of  p  which  are  off-diagonal  in  the  island  charge 
index  n—  there  is  no  coherence  between  different  charge  states,  as  tunneling 
events  necessarily  create  an  electron-hole  excitation.  Further,  if  we  focus 
on  small  qubit  frequencies,  we  may  continue  to  restrict  attention  to  only 
n  =  0  and  n  =  1  (i.e.  H  is  not  large  enough  to  “turn  on”  tunneling 
processes  which  are  normally  energetically  forbidden).  Thus,  there  are  8 
relevant  matrix  elements  of  p—  for  each  of  the  four  qubit  density  matrix 
elements  (i.e.  gg,  ee,  ge,  eg),  there  are  two  possible  island  charge  states. 
We  combine  these  elements  into  a  vector  p  =  (pgg,  pee,  pge,  peg),  where 

Pag  =  (<0,5o|/9|0,5o),<l,5i|p|l,5i>),etc.  Organizing  the  resulting  evolution 
equation  in  powers  of  the  coupling  A,  we  obtain  in  the  Schrodinger  picture: 

j/=  (Ao  +  ^Ai  +  ^A2  +  ...)p  (63) 
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We  discuss  the  significance  of  the  matrices  Aj  in  what  follows. 

The  8x8  matrix  Ao  describes  the  evolution  of  the  system  at  zero 
coupling: 


Aq 


/MO  0  0  \ 

0  M  0  0 

0  0  +iQ  +  M'  0 

\  0  0  0  -in  +  M'  ) 


(64) 


with  the  2x2  matrices  M  and  M'  being  defined  by: 


M  = 


-(rlo  +  r/o) 

rib  +  rio 


roi  +  roi  A 
-(r0i  +roi)  ) 


(65) 


At  zero  coupling  there  are  no  transitions  between  different  qubit  states,  and 
hence  Ao  has  a  block-diagonal  form.  There  are  two  independent  stationary 
solutions  of  Eq.  (63)  at  A  =  0  (i.e.  two  zero  eigenvectors  of  Ao),  which 
correspond  to  being  either  in  the  qubit  ground  or  qubit  excited  state: 


Zg  =  (p0, PI, 0,0, 0,0, 0,0)  , 


ze  =  (0,0,po,pi,0,0,0,0)  . 


(66) 


(Po,Pi)  are  the  stationary  probabilities  of  being  in  the  n  =  0  or  n  =  1 
charge  states: 


bo,  Pi) 


01 


10 


+  Tio  ’  Tqi  +  r 


10 


(67) 


The  existence  of  two  zero-modes  is  directly  related  to  the  fact  that  at  zero 
coupling  (A  =  0),  the  probabilities  to  be  in  the  qubit  ground  and  excited 
state  are  individually  conserved. 

At  non-zero  coupling,  the  matrices  Ai  and  A2  appearing  in  Eq.  (63) 
generate  transitions  between  different  qubit  states.  The  matrix  A2  directly 
couples  pgg  and  pee,  while  Ai  couples  pgg  and  pee  to  the  off-diagonal  blocks 
pge  and  peg.  The  effect  of  these  matrices  will  be  to  break  the  degeneracy  of 
the  two  zero  modes  of  Eq.  (63)  existing  at  A  =  0.  After  this  degeneracy  is 
broken,  there  will  still  be  one  zero  mode  po,  describing  the  stationary  state 
of  the  coupled  system  (the  existence  of  a  stationary  solution  is  guaranteed 
by  the  conservation  of  probability).  For  weak  coupling,  the  qubit  density 
matrix  obtained  from  po  will  be  diagonal  in  the  basis  {|p(n)),  |e/n\)},  which 
corresponds  to  the  average  SET  charge  (n)  =  p\ .  The  ratio  of  the  occu¬ 
pancies  of  these  two  qubit  states  will  yield  the  ratio  between  the  relaxation 
rate  Tj_  and  the  excitation  rate  Tt.  In  addition,  there  will  also  be  a  slow, 
time-dependent  mode  of  Eq.  (63)  arising  from  breaking  the  degeneracy  of 
the  two  A  =  0  zero  modes.  This  time-dependent  mode  will  describe  how 
a  linear  combination  of  zg  and  ze  relaxes  to  the  true  stationary  state,  and 
will  have  an  eigenvalue  A  =  — Tt  —  Ti ,  i.e.  the  mixing  rate. 
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Thus,  we  need  to  do  degenerate  second  order  perturbation  theory  in 
the  coupling  A  to  obtain  the  relaxation  and  excitation  rates  Tj  and  Tp 
The  only  subtlety  here  is  that  the  matrix  M  is  not  Hermitian,  implying 
that  it  has  distinct  right  and  left  eigenvectors.  Letting  z  represent  the  left 
eigenvector  of  Ao  corresponding  to  the  right  eigenvector  z,  we  define  the 
projector  matrix  P  as: 

P  =  \zg){*g\  +  \ze)(Ze\,  (68) 

and  let  Pj_  denote  1  —  P.  As  usual,  degenerate  second  order  perturba¬ 
tion  theory  requires  diagonalizing  the  perturbation  in  the  space  of  the 
degenerate  eigenvectors.  We  are  thus  led  to  look  at  the  matrix  Q,  defined 
as: 

Q  =  ^(PA2P  +  PA1P±[-A0]-1P±A1P)  (69) 

/'.From  the  definition  of  Q,  we  may  immediately  identify  the  rates  Ft  and 
ri: 

rT  =  (ze\Q\zg)  Tj  =  (zg\Q\ze)  (70) 

We  thus  see  how  the  rates  Tt  i  arise  in  the  present  approach-  they 
are  related  to  breaking  the  degeneracy  between  two  zero-modes  (stationary 
solutions)  which  exist  at  zero  coupling.  Note  that  there  are  two  distinct 
contributions  to  Tt  i,  coming  from  the  two  terms  in  the  matrix  Q:  a  “direct” 
contribution  involving  A2  and  an  “interference”  contribution  involving  Ai 
acting  twice.  These  two  terms  have  a  different  physical  interpretation,  as 
will  become  clear. 

Let  us  first  consider  the  rate  Tt,  which  describes  how  noise  in  the  SET 
causes  ground  to  excited  state  transitions  in  the  qubit.  For  this  rate,  our 
approximation  of  only  keeping  two  charge  states  will  be  valid  for  all  splitting 
frequencies  12.  To  evaluate  the  “direct”  contribution  to  this  rate,  which 
involves  the  first  term  in  the  matrix  Q,  note  that  the  relevant  part  of  A2 
has  the  expected  form: 

=  (  r°f0  ro"  )  <71» 

i.e.  it  consists  of  tunnel  rates  which  correspond  to  having  given  up  an  energy 
12  to  the  qubit.  Using  Eqs.  (69)  and  (70),  we  find: 

rrl—  =  (£)  Vr„+«rM)  =  (A)2  (EE£h±Ii2£h)  (72) 

The  direct  contribution  to  Tt  has  a  very  simple  form:  for  each  charge  state 
n  =  0,1,  add  the  rate  to  tunnel  out  of  n  while  exciting  the  qubit,  weighted 


rjsdraftl2.tex;  11/10/2002;  0:40;  p.25 


26 


by  both  the  probability  to  be  in  state  n,  and  the  overlap  between  ground 
and  excited  states  (i.e.  (A/fi)2).  This  is  very  similar  to  how  one  typically 
calculates  the  current  for  a  SET:  one  adds  up  the  current  associated  with 
each  charge  state  (i.e.  a  difference  of  rates),  weighted  by  the  occupancy 
of  the  state.  The  direct  contribution  to  Tj  neglects  any  possible  coherence 
between  successive  excitation  events;  as  a  result,  it  fails  to  recover  the 
classical  expression  of  Eq.  (54)  in  the  small-O  limit. 

We  now  consider  the  “interference”  contribution  to  Tj  coming  from 
the  second  term  in  the  expression  for  matrix  Q  (c.f.  Eq.  (69)).  After  some 
algebra,  we  obtain  the  following  for  the  interference  contribution  to  Tt: 


2A2 

Ttlmt  =  ~~qT  (^r10  +  Tir01) 


(rs)2 

n2  +  (ts)2 


(73) 


where: 


rs  = 


Fro  +  r+  +  r, 


10 


01 


(74) 


This  contribution  is  purely  negative,  and  is  only  significant  (relative  to 
the  direct  contribution)  at  low  frequencies  Q  <  T.  We  can  interpret  this 
equation  as  describing  the  interference  between  two  consecutive  excitation 
events.  For  example,  consider  the  first  term  in  Eq.  (73).  This  describes 
a  process  where  a  SET  initially  in  the  charge  state  n  =  0  undergoes  a 
tunnel  event  to  the  n  =  1  state,  creating  a  superposition  of  qubit  ground 
and  excited  states.  At  some  later  time  the  SET  relaxes  to  the  station¬ 
ary  distribution  (po,pi)  of  the  charge  states,  again  partially  exciting  the 
qubit.  Letting  At  represent  the  time  between  these  two  events,  we  have  the 
approximate  sequence: 

|0, go)  ^  |l,5i>  +a|l,ei) 

— — >  emAt/2\l,gi)+e-inAt/2a\l,ei) 

(eiClAt/2/3  —  e~inAt/2a)  |0,  e0)  + 


(75) 

(76) 


Here,  a  is  the  amplitude  associated  with  qubit  excitation  having  occurred 
during  the  first  (n  =  0  — >  1)  tunnel  event,  while  (3  is  the  amplitude 
associated  with  excitation  occurring  during  the  second  (n  =  1  — >  0)  tun¬ 
neling.  These  amplitudes  will  be  given  by  the  corresponding  matrix  overlap 
elements: 

a  =  (eibo)-^  P=(e0 \gi)  ~ (77) 

In  the  final  state  after  the  two  tunnelings  (Eq.  (76)),  there  are  two  terms 
in  the  amplitude  of  the  state  |0,  eo),  corresponding  to  the  fact  that  qubit 
excitation  could  have  occurred  in  either  the  first  or  the  second  tunnel  event. 
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To  get  a  rate  for  this  double  excitation  event,  we  should  take  the  modulus 
squared  of  the  final  1 0,  eo)  state  amplitude,  then  multiply  by  the  occupancy 
of  the  initial  state  (po)  and  the  rate  of  the  first  tunnel  event  (Tio).  The 
interference  term  in  the  resulting  expression  takes  the  form: 

2A2 

rT|mt  =  (porio)  x  2Re  ( a*(3elQAt )  =  -  (p0T w) cos(Q, At)  (78) 

The  above  expression  is  a  function  of  the  time  At  between  the  first  and  sec¬ 
ond  tunnel  events.  This  time  is  determined  by  the  fact  that  the  intermediate 
superposition  state  (Eq.  (75))  corresponds  to  a  non-stationary  distribution 
of  charge  on  the  SET  island,  and  will  decay  via  tunneling  to  the  stationary 
distribution  (po,Pi)  at  a  rate  Ts-  Taking  this  decay  to  be  Poissonian,  and 
averaging  over  At,  we  obtain: 


r-rlmt  —  — (pobio)- 


n2  n2  +  (Ts)2 


(79) 


This  is  precisely  the  first  term  in  Eq.  (73);  the  second  term  can  be  ob¬ 
tained  in  the  same  way,  by  now  considering  a  situation  where  the  SET  is 
initially  in  the  n  =  1  charge  state.  As  claimed,  T y | corresponds  to  the 
interference  between  two  consecutive  excitation  events.  The  negative  sign 
of  this  contribution  can  be  directly  traced  to  the  matrix  overlap  elements 
(c.f.  Eq.  77).  Also,  we  see  that  the  suppression  of  the  interference  term  at 
large  hi  results  from  phase  randomization  occurring  during  the  delay  time 
between  the  two  excitation  events. 

Returning  to  the  total  noise,  we  combine  Eq.  (73)  with  the  direct  con¬ 
tribution  Eq.  (72)  to  TT;  comparing  against  Eq.  (57),  we  obtain  the  final 
expression  for  Sq(Q.)  at  all  negative  frequencies: 


Sq(-M)  = 


P0ri0  +PlTn 


+  I  (rio  +  rio  +  roi  +  roi)" 


Note  for  large  |fl|  (i.e.  |fl|  >  max(A£'g1,  AE“0)  ~  Vds/ 2),  Sq(—Q)  will 
vanish  identically  at  zero  temperature.  Physically,  this  cutoff  corresponds 
to  the  largest  amount  of  energy  the  SET  can  give  up  to  the  qubit  during  a 
single  tunnel  event;  giving  up  more  energy  would  suppress  the  event  com¬ 
pletely  (i.e.  the  tunnel  rates  have  a  step-function  form  at  zero  temperature, 
c.f.  Eq.  (50)).  If  one  included  higher  order  processes  in  the  tunneling  (i.e. 
went  beyond  sequential  tunneling),  correlated  tunneling  events  involving 
the  full  voltage  drop  over  both  junctions,  Vds-,  would  move  this  cutoff  to 
higher  values  of  absolute  frequency. 

We  now  turn  to  the  relaxation  rate  T  i ,  and  hence  the  positive  frequency 
parts  of  Sq.  The  calculation  proceeds  exactly  as  that  for  r p ,  the  only 
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Figure  5.  Symmetrized  SET  charge  noise  as  a  function  of  frequency,  for  typical  SET 
parameters:  g  =  1,  Ec/Icb  =  2K,  J\f  =  0.33,  cVds  =  Ec,  and  T  =  20m K.  The  dashed 
line  is  the  classical  telegraph  noise  (Eq.  (54)),  while  the  dot-dashed  line  is  the  noise  of 
two  parallel  tunnel  junctions  (Eq.  (55)).  Inset:  full  (non-symmetrized)  quantum  noise  for 
identical  SET  parameters;  the  dashed  line  is  the  symmetric  classical  telegraph  noise. 


modification  being  that  one  now  needs  to  include  the  charge  states  n  =  2 
and  n  =  —1,  as  the  SET  could  absorb  enough  energy  from  the  qubit  to 
make  transitions  to  these  states  possible.  We  can  combine  the  result  for  Tj_ 
with  Eq.  (80)  to  obtain  a  single,  compact  expression  for  the  noise  at  all 
frequencies  first  obtained  by  Johansson  et.  al  [14]:  13 

s  ,  n  =  pq  FioH  +  r_i,o(q;)]  +pi  [TqiH  +  TgiM]  . 

u2  +  \  [rio(w)  +  Tio(— ui)  +  Toi(w)  +  r0i(-w)]2 

Shown  in  Figure  5  is  the  symmetrized  noise  Sq(uj)  +  Sq(-lu)  for  typical 
SET  parameters.  One  can  clearly  see  abrupt  changes  in  the  slope  of  this 
curve;  each  of  these  kinks  corresponds  to  a  threshold  frequency  at  which  a 
given  tunneling  process  either  turns  on  or  turns  off.  For  comparison,  curves 

13  Eq.  (81)  ignores  additional  order  g/( 2n)  terms  which  arise  in  the  denominator  at 
positive  frequencies  large  enough  to  turn  on  tunneling  to  higher  charge  states;  such  terms 
are  clearly  negligible  in  the  sequential  tunneling  regime  due  to  the  smallness  of  g/(2n). 
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Figure  6.  a)  Average  charge  state  of  a  Cooper  pair  box  coupled  to  a  SET,  as  a  function 
of  box  gate  voltage,  using  identical  SET  parameters  as  above.  The  box  parameters 
are  Ec/kB  =  0.5K,  Ej/Ub  =  0.25K  and  the  coupling  constant  is  k  =  0.04.  We  also 
include  relaxation  effects  due  to  a  10%  coupling  to  a  50S7  environment.  The  dashed  curve 
corresponds  to  assuming  the  SET  produces  classical  telegraph  noise,  the  solid  curve 
corresponds  to  using  the  full  quantum  noise  of  the  SET,  and  the  dashed-dot  curve  is  the 
box  ground  state,  b)  The  relaxation  time  Ti  for  the  same  system,  as  a  function  of  box 
gate  voltage. 


corresponding  to  classical  telegraph  noise  and  to  the  uncorrelated  noise  of 
two  tunnel  junctions  are  also  shown.  At  low  frequencies  the  symmetrized 
true  noise  matches  the  classical  curve;  for  higher  frequencies,  it  lies  above 
the  classical  curve  but  below  the  curve  corresponding  to  the  uncorrelated 
case.  The  inset  of  this  figure  shows  the  both  the  negative  and  positive 
frequency  parts  of  Sq(u). 

It  is  easy  to  check  that  in  the  limit  uo  — >  0,  Eq.  (81)  recovers  the  classical 
telegraph  expression  of  Eq.  (54).  In  the  high-frequency,  zero  temperature 
limit,  one  can  also  see  that  Eq.  (81)  approaches  the  uncorrelated  result  of 
Eq.  (55)  from  below: 


Sq(u) 


©M— 

7 TCO 


(82) 


Note  that  at  high  frequencies,  it  is  only  the  “direct”  terms  which  contribute 
to  the  noise-  the  interference  contribution  is  not  important  in  the  limit 
of  uncorrelated  tunneling.  The  fact  that  the  noise  approaches  the  high 
frequency  limit  from  below  results  from  the  tendency  of  charging  energy 
induced  correlations  (which  are  present  for  a  finite  uj/Eq)  to  suppress 
fluctuations  of  n,  and  thus  suppress  the  noise.  Note  that  the  interpolation 
between  the  low  and  high  frequency  limits  here  is  very  different  than,  e.g., 
interpolating  between  thermal  noise  and  zero-point  fluctuation  noise  in  a 
tunnel  junction.  In  the  latter  case,  one  is  effectively  combining  two  sources 
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of  noise;  here,  one  is  simply  turning  off  correlations  brought  on  by  the 
charging  energy  by  increasing  lo. 

Finally,  shown  in  Figure  6  is  the  average  charge  state  of  a  Cooper  pair 
box  coupled  to  a  SET  with  identical  parameters  to  that  in  Fig.  5.  We  have 
also  included  here  the  relaxation  effects  of  the  environment,  modelled  as 
in  Sec.  6  as  a  50fl  impedance.  Note  that  even  near  the  box  degeneracy 
point,  there  are  large  deviations  between  the  result  obtained  using  the  full 
quantum  noise  of  the  SET  and  that  obtained  from  using  only  classical 
telegraph  noise.  In  Fig.  6b,  we  show  the  relaxation  rate  T\  for  the  same 
system.  Note  that  the  differences  between  using  the  full  quantum  noise  and 
the  classical  expression  are  not  so  evident  here. 

9.  Summary 

In  this  article,  we  have  emphasized  the  need  to  discuss  quantum  noise 
processes  using  their  two-sided  spectral  densities.  Because  of  the  quantum 
nature  of  noise,  the  positive  and  negative  frequencies  are  generally  unequal, 
in  order  to  account  for  spontaneous  emission.  A  two-level  system  was  shown 
to  be  an  ideal  spectrum  analyzer  for  probing  the  quantum  nature  of  a 
noise  process  or  reservoir.  With  the  advent  of  real  electrical  circuits  which 
behave  as  coherent  two- level  systems  (e.g.,  [5], [7]),  we  can  now  build  and 
use  quantum  electrical  spectrum  analyzers.  We  also  described  the  use  of 
a  qubit  as  a  theoretical  tool,  by  following  the  evolution  of  the  density 
matrix  of  a  TLS  coupled  to  the  noise-producing  system  of  interest.  This 
technique  appears  to  be  quite  powerful,  as  it  can  yield  analytical  results 
for  the  full  quantum  noise  spectrum  of  a  wide  variety  of  devices,  including 
the  superconducting  SET  [17].  The  distinction  between  the  classical  noise 
and  the  quantum  noise,  found  in  this  way,  leads  to  dramatically  different 
predictions  (c.f.  Fig.  6  and  Ref.  [17])  for  continuous  measurements  of  qubits 
with  an  SET.  The  “coupled-system”  calculational  approach  also  allows 
predictions  of  the  dephasing  by  the  measurement,  the  performance  relative 
to  the  Heisenberg  uncertainty  limit  [13],  the  fidelity  of  single-shot  measure¬ 
ments  of  the  qubit  states,  and  the  effects  of  strong  coupling  to  the  qubit. 
The  combined  theoretical  and  experimental  advances  raise  many  interesting 
possibilities  for  testing  our  understanding  of  quantum  measurement  theory 
with  mesoscopic  devices. 
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